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This settles positively the 3-dimensional Propeller Conjecture of Khot and Naor (FOCS 2008). The 
proof of m reduces the problem to a finite set of numerical inequalities which are then verified 
with full rigor in a computer-assisted fashion. The main consequence (and motivation) of ([T]) is 
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4x4 centered and spherical hypothesis matrix equals 
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Figure 1. The partition o/M'^ that maximizes the sum of squared lengths of Gauss- 
ian moments is a "propeller": three planar 120° sectors multiplied by an orthogonal 
line, with the rest of the partition elements being empty. 
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1. Introduction 



For X = {xi, . . . ,Xm),y 
product, and let ||x||2 = \/ {x, x) 
X = {xi,X2, X'i),y = (yi, 2/2,2/3) G 
Gaussian measure on I 



(yi, . . . , Um) G K"* let {x, y) = YllLi ^iyi denote their standard scalar 
denote the corresponding Euclidean norm. The cross product of 
is denoted xxy = {x2yz-x^y2,xzyi-xiy-i,xiy2-X2yi). The 
i.e., the measure whose density is x 1— )• (27r)^'"/^e~ll^ll2/^, is denoted 7™,- 



The following theorem is our main result, asserting that among all measurable partitions of M , 
the "propeller partition" as depicted in Figure [T] maximizes the sum of the squared lengths of the 
Gaussian moments associated to the members of the partition. 

Theorem 1.1 (Main theorem; geometric formulation). Let {^1, . . . , ^fc} be a partition o/M'^ into 
Lebesgue measurable sets. For i G {1, . . . , A;} let Q = X4. xd'y^lx) G be the Gaussian moment of 
the set Ai. Then 

k 



(2) 



cannot 



Let {Pi, P2, P3} be the partition o/M^ into 120° sectors centered at the origin. The bound 
be improved, with equality holding if Ai = Pi for i G {1, 2, 3} and Ai = 9 for i £ {4, . . . , k} (up to 
measure zero corrections, orthogonal transformations, and renumbering of the sets {Ai, . . . , Aj^} ). 

The Propeller Conjecture was posed by Khot and Naor in [23] as part of their investigation of the 
computational complexity of the Kernel Clustering problem from Machine Learning. Specifically, 
they conjectured the validity of the bound ^ for measurable partitions of M"^ for all m > 2. They 
proved this conjecture for m = 2, and for m = 1 they showed that the right hand side of ([2]) can be 
improved to the sharp bound ^. It is also shown in |23j that proving ([2]) for k = A implies the same 
conclusion for all A; G N. Correspondingly, for partitions of M"^ it suffices to prove this statement 
for k = m + 1. 

Before explaining the complexity-theoretic consequence of Theorem 11.11 which was the motiva- 
tion of [23] for posing the Propeller Conjecture, we state two equivalent formulations of it: the first 
probabilistic (a sharp estimate for the expected maximum of a Gaussian vector) and the second 
analytic (a sharp Grothendieck inequality). The equivalence of these results to Theorem 11.11 was 
established in [23]; for Theorem 11.21 below see [23^ Lem. 3.7] and for Theorem 11.31 below see [23^ 
Thm. 1.1]. 



Theorem 1.2 (Main theorem; probabilistic formulation). Let {gi, §2, 93, 94,) G 
Gaussian vector (with arbitrary covariance matrix), then 
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max gi 
ie{l,2,3,4} 



4 



be a mean zero 



(3) 



The bound ([3]) cannot be improved, with equality holding when the covariance matrix of {gi, g2, 93, g^) 
equals 

/ 2 -1-1 0\ 



-1 








-1 2 
-1 -1 
\ 

The following analytic formulation of Theorem 1 1.1 1 is in terms of a sharp Grothendieck inequality. 
Grothendieck [13] pioneered in 1953 the use of inequalities of this type in functional analysis, and 
ever since then such inequalities have permeated many mathematical disciplines. See [23, 24, ,25] for 
an explanation of how the result below relates to a natural extension of the Grothendieck inequality 
and to approximation algorithms. 
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Theorem 1.3 (Main theorem; analytic formulation). Let [aij) he an n ^ n positive semidefinite 
matrix with '^^^i Sj=i '^u ~ ^- ^^"^ every {^1,^2,^3, ^4} ^ S'^ with Ylt=i Vi = we have 



max V V aij{xi,Xj) — max V V aij{yi, yj). 

3;i,...,x„e5"-i ■'— ■ ■'^ — ; o »/i,.--,j/nG{i)i,«;2,'y3,^'4} — : — : 

2=1 J = l 2=1 7 = 1 



(4) 



The bound (jl]) cannot be improved, with asymptotic equality holding when {fi, ^2, ^"3, ^4} are the 
rows of the following matrix. 
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5*66 [241 Sec. 3.1] for a description of a family of n x n matrices (ajj) for which equality in @ is 
asymptotically attained (as n ^ 00). 

Inequality (jl]) was proved in [23]. Our new contribution is the assertion that @ cannot be 
improved, a statement that is equivalent to inequality ^ of Theorem 1 1.1[ This is not the first time 
that attempts to prove sharpness of a Grothendieck inequality led to an extremal geometric par- 
titioning question. Notably, see Konig's conjecture [27] as a step towards Krivine's conjecture [29] 
on the sharpness of his version of the classical Grothendieck inequality. Unlike the Propeller Con- 
jecture in M^, these conjectures turned out to be false [8], but they do indicate the interconnection 
between Grothendieck inequalities and extremal geometric partitioning problems. 

The main consequence (and motivation) of Theorem II. II is complexity theoretic. To explain it we 
briefly recall Khot's Unique Games Conjecture (UGC), which asserts that for every e S (0, 1) there 
exists a prime p = p{e) S N such that no polynomial time algorithm can perform the following 
task. The input is a system of m linear equations in n variables xi, . . . ,Xn, each of which has 
the form Xi — Xj = Cij mod p. The algorithm must determine whether there exists an assignment 
of an integer value to each variable Xi such that at least (1 — e)m of the equations are satisfied, 
or no assignment of such values can satisfy more than em of the equations. If neither of these 
possibilities occurs, then an arbitrary output of the algorithm is allowed. The UGC was introduced 
by Khot in [20], though the above formulation of it, which is equivalent to the original one, is 
due to [22]. The use of the UGC as a hardness hypothesis has become popular over the past 
decade; we refer to Khot's survey [21] for more information on this topic. Saying that the UGC 
hardness threshold of an optimization problem equals a € [l,c>o) means that for every e £ (0, 1) 
there exists a polynomial time algorithm that outputs a number that is guaranteed to be within a 
multiplicative factor of a -|- e from the solution of ff, and that the existence of such an algorithm 
with approximation guarantee of a — e would contradict the UGC. 

The Kernel Clustering problem is a clustering framework for covariance matrices that originated 
in the work of Borgwardt, Gretton, Smola and Song ^36j in the context of Machine Learning. This 
problem is generic in the sense that it contains well-studied optimization problems as special cases, 
and its versatility allows one to design a variety of algorithms tailor-made for particular applications 
(many of these algorithms are at present shown [36] to be successful empirically, but not rigorously). 
The input of the Kernel Clustering problem is an n x n symmetric positive semidefinite matrix 
A = (oij) and a. k x k symmetric positive semidefinite matrix B = (bij) called the hypothesis 
matrix. Think of n as very large and k as small, the goal being to cluster the entries of A into 
a k X k matrix that is most correlated with the hypothesis matrix B. Formally, given a partition 
{Si, . . . , S"/;} of {1, . . . , n}, form the associated clustered version of A by summing the entries of A 
over the blocks induced by the partition {Si, . . . , Sk}- One thus obtains a k x k matrix C = (cij) 
given by Cij = Y^^^ t)(^SixSj ^^t- Let Clust(j4|i?) denote the maximum of J2i=i Sj=i '^ij^ij 
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partitions {Si, . . . , Sk} of {1, ... , n}. We refer to [361 ESI 1211 ES] for further explanation of this 
clustering framework, as well as a discussion of important special cases arising from appropriate 
choices of the hypothesis matrix B. 

In what follows, an n x n symmetric positive semidefinite matrix A = {aij) is called centered if 
Yll=i X]j=i ^ij ~ ^^'^ called spherical if an = 1 for alH G {1, . . . , n}. 

Theorem 1.4 (Main theorem; complexity theoretic formulation). Let G he the following optimiza- 
tion problem. The input is an nx n symmetric positive semidefinite centered matrix A = (aij), and 
also a 4 X 4 symmetric centered and spherical positive semidefinite matrix B = [hij). The goal is 
to compute in polynomial time the quantity C\\ist{A\B). Then the UGC hardness threshold of G 
equals 

A polynomial time algorithm with approximation ratio ^+o(l) for the problem of Theorem ll.4l 
was designed in [23] . The fact that Theorem 11.11 implies the matching UGC hardness result was 
also proved in [23]. More generally, it was shown in [23] that the validity of the m-dimensional 
Propeller Conjecture for some m ^ 3 would imply that the UGC hardness threshold of the variant 
of G with B being an m x m matrix equals ^ (l — 

This is not the first time that extremal problems in the measure space (M'",7m) arose from 
investigations in complexity theory; see for example the Majority is Stablest Conjecture of Khot, 
Kindler, Mossel and O'Donnell [2^, which was solved by Mossel, O'Donnell, and Oleszkiewicz |33] 
via a reduction to a classical isoperimetric inequality of Borell [7] . A special feature of the Propeller 
Conjecture is that its conclusion is quite surprising: despite allowing for partitions of into m+1 
sets, the optimal partition has only three nonempty sets. Note that this degeneracy property occurs 
for the first time when m = 3, since in the result of [23] for partitions of into three sets, all 
three sets are indeed present in the optimal partition. Thus Theorem 11.11 is a proof of the first 
nonintuitive case of the Propeller Conjecture: a verification of a prediction about a clean yet 
unexpected geometric phenomenon that arose from an investigation in computational complexity. 
To the best of our knowledge this is the first time that such a development occurred in this theory. 

Our proof of Theorem 11.11 proceeds as follows. First, using reductions of the Propeller Conjecture 
that were obtained in j23j . we show that it suffices to prove an analogous partitioning problem for 
the sphere S'^. Namely, given a partition of into four spherical triangles, the goal is to show that 
the sum of the squares of lengths of the moments (with respect to the surface area measure on S^) 
of these spherical triangles cannot exceed 97r^/4. Using additional geometric arguments, we further 
reduce the question to an optimization problem over a three dimensional search space. We thus 
obtain a coordinate system with three degrees of freedom, with respect to which one must solve 
a rather complicated nonlinear optimization problem. After proving several additional estimates 
that serve to further reduce the search space and give crucial modulus of continuity estimates, 
we show that it suffices to check that the desired estimate holds true for a finite list of spherical 
partitions (arising from an appropriate net of the search space) . This list of inequalities that must 
be proved is explicit but very large, so we proceed to check it in a fully rigorous computer-assisted 
fashion. That is, our computation carefully accounts for all the rounding errors; equivalently the 
computation can be viewed as an implementation in our setting of interval arithmetic (see |17|). 

The role and implications of rigorous computer-assisted proofs has been discussed at length in 
the literature; many such discussions appear in papers that establish striking results via a proof 
that has a computer-assisted component (famous examples include [H [21 [181 [TTl [151 [12] ) . We see 
no reason to include here a new treatment of this topic. Zwick's discussion in |39j is an excellent 
reference for a well thought out explanation of the role of computer-assisted proofs in mathematics 
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and computer scienc^H- The essence of the argument can be best conveyed by quoting Zwick 
directly |39l Sec. 7]: 

"A typical computer assisted proof, like ours, is essentially composed of two steps. 
The first step says something like: "Here is program P. If program P, when executed 
by an abstract computer, outputs YES, then the theorem is true, because . . . ". 
The second step says something like: "I ran program P and it said YES!", or more 
specifically: "I compiled program P using compiler A, under operating system B, 
and ran it on processor C, again under operating system B. It said YES!"." 

Zwick then proceeds to explain that 

"The first step is a conventional mathematical proof. It simply argues that a certain 
program has certain properties. Such arguments are common in computer science. 
Anyone who objects to this step of the proof should find a flaw, or a gap, in the 
arguments made." 

Nevertheless, Zwick explains that the second step, i.e., the claim that the compilation of the program 
gave the desired result, is not of the same nature: 

"The second step is more problematic. It is certainly not a proof in the conventional 
mathematical sense of the word. Many things can go wrong here. Program P may 
not have been compiled correctly by compiler A. Or, due to some bug in operating 
system B, the program did not run as intended. Or, processor C may suffer from a 
design flaw that causes an incorrect execution of the program, or perhaps the specific 
processor used has a short-circuit somewhere, and so on. It is possible to reduce 
the likelihood of such problems by compiling the program using several different 
compilers, and running it on different processors. But, while we may eventually be 
able to produce mathematical correctness proofs for compilers, operating systems, 
and the hardware design of the processors, it seems that we would never be able 
to produce a mathematical proof that no hardware fault occurred during a specific 
computation. 

What we can do, is instruct the computer executing program P to print a trace 
of the execution. This trace is a mathematical proof, though perhaps a not so 
inspiring one. Like any mathematical proof it should be checked carefully to make 
sure that it is correct. In many cases, however, this is not humanly possible. The 
main problem with computer assisted proofs, therefore, is that they are usually too 
long. They are also, in most cases, less insightful than conventional proofs." 

The bulk of the work presented in this article consists of conventional geometric and analytic 
proofs, resulting in a reduction of the problem to a region where the desired inequality holds with 
"room to spare". This extra room allows us to complete the proof by checking a finite list of 
concrete inequalities (that we write explicitly) on a certain finite net. Each such inequality can be 
checked by hand, but the number of such checks is large, and it would be unrealistic (and probably 
unilluminating) to complete this final check manually. The code that we used is publicly available 
at the following link, which contains a command-line user interface so that interacting with the 
code is easy. 

http : //cims . nyu . edu/~aukosh/propeller/propeller . html 



Zwick's paper [39] is another instance of a computer-assisted proof in the context of approximation algorithms. 
While the present article deals with an integrality gap lower bound for a semindefinite program (the left hand side 
of (|4])), Zwick proves an integrality gap upper bound, i.e., he uses a rigorous computer-assisted argument to show 
that a certain algorithm performs well rather than to prove a hardness result. 
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The code for this proof was run many times on several different computers. The technical require- 
ments are quite modest, with the processor speed only affecting the total run time, and the total 
memory requirement being less than 200 Mb. One can run the code on anything from a dedicated 
compute server to a laptop computer (we did both). We ran the code in two different ways. The 
simplest way was to run it in a serial fashion. That is, we ran the program as a single instance 
which performed each step of the algorithm described in this paper one after the other. This took 
37 hours to complete on a single processor of a compute server with 3.3 GhZ Intel Xenon processors. 
The second way we ran it was to parallelize the algorithm by breaking up the domain of interest 
into seven distinct portions and running separate instances on each. This procedure took eight 
hours when split between seven 3.3 GhZ Intel Xenon processors. (We also ran the program split 
between seven 2.4 GhZ AMD Opteron processors; this procedure took twelve hours.) 

Needless to say, our proof of the Propeller Conjecture in leaves something to be desired, 
since we do not have a short explanation of the validity of the computations in the final step of 
our proof. It is conceptually important to have a proof of this result, despite the fact that it 
ends in a lengthy computation, not only because it yields interesting results in mathematics and 
computational complexity. Importantly, since the Propeller Conjecture is a nonintuitive prediction 
that arose from investigations in computational complexity, knowing that the first nonintuitive 
case of this conjecture is indeed true puts the full Propeller Conjecture, and consequently the link 
between geometry and algorithms that it describes, on a stronger footing. Thus, in addition to 
yielding a remarkably involved UGC hardness result, our theorem will hopefully invigorate future 
research on this problem that might lead to a traditional mathematical proof of the Propeller 
Conjecture in that can hopefully be extended to higher dimensions. 

A key feature of the Propeller Conjecture is that it asserts that a natural optimization problem 
exhibits intermediate-dimensional symmetry breaking. There are precedents of results of this type 
that have been proved in the literature, mostly using Fourier-analytic methods. For this reason we 
are hopeful that a clean and shorter proof of the Propeller Conjecture will eventually be found. 
Consider for example the problem of sharp Khinchine inequalities: if ei, . . . , e„ are i.i.d. symmetric 
random variables taking values in { — 1, 1} and p € [1,2), the goal is to compute the minimum of 
E ajEjl^] over all unit vectors a = (ai, . . . , a„) G M". For p = 1 it was conjectured by Little- 

wood (see [16]) that this minimum occurs at a = (1, 1, 0, ... , 0)/\/2. Littlewood's conjecture was 
solved affirmatively by Szarek [37] (see also [381 EO] ) . Haagerup [TJ] proved that the same two di- 
mensional symmetry breaking occurs for p S [IjPo]) where po = 1.87... is the solution of the equation 
2r((p + l)/2) = ^/^T, i.e., the unit vector that minimizes E [|^"^-^ ajej|^] equals (1,1,0,... ,0)/\/2 
for p G [1,Pq] and for p G [po, 2] this minimum occurs at a = (1, ... , l)/y^. Another famous result 
of this type is Ball's cube slicing theorem [3] (resolving a conjecture of Hensley [H]), asserting 
that the hyperplane section of [—1, 1]"' with maximal (n — l)-dimensional volume is perpendicular 
to (1, 1, 0, . . . , 0); see [35] for the corresponding result for complex scalars, as well as [H E] for 
an analogous statement for projections. We refer to [34] for a unified treatment of the results of 
Szarek, Haagerup, and Ball. There is also a conjecture of Milman predicting a symmetry breaking 
phenomenon for extremal volumes of slabs in the cube [—1, 1]"; see [5l[28] for partial results along 
these lines. We mentioned the above statements since we believe that among the literature they 
are most similar to the Propeller Conjecture, and perhaps (with much more work) related methods 
could be be used to address the Propeller Conjecture as well. 

This article is organized as follows. In Section [2] we present the reduction of the Propeller 
Conjecture to a certain optimization problem for spherical partitions, and we explain the main 
ingredients of our proof, including the conventional geometric/analytic arguments, as well as the 
computer-assisted component. The proofs of the geometric and analytic results leading to the 
final computational step are contained in Section [3] and Section HI Section 14.11 contains a detailed 
explanation of how the numerical step is implemented so as to account for all possible rounding 
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errors. We remark that our arguments extend mutatis mutandis to higher dimensions. For the sake 
of simplicity we present the entire argument in M^, since while it is conceivable that the scheme 
presented here can yield a computer-assisted proof of the Propeller Conjecture in higher dimensions, 
at some fixed dimension the computer-assisted component of the proof will become unfeasible. Now 
that we know that a nonintuitive case of the Propeller Conjecture is indeed correct, the next natural 
step is to search for a proof that extends to all dimensions, rather than attempting to prove a few 
more cases in low dimensions. 



2. An overview of the proof of Theorem 11.11 

From now on assume for the sake of eventually obtaining a contradiction that {Ai}j^^ is a 
partition of M'^ into measurable sets that violates the propeller conjecture, i.e.. 



4 

E 

i=l 



xd-f3{x) 
A, 



9 

> i- (5) 

Sir 



Assume moreover that the maximum of X]i=i Jb ^^73 (^) over all measurable partitions {Bi}f^^ 

of is attained at {Ai}j^i. For a proof that this maximum is indeed attained, see \23\ Lem. 3,1]. 
Using this maximality, it follows from \23\ Lem. 3.3] that (up to measure zero corrections) the Ai 
are cones with cusp at the origin, and if we write d = xd'y3{x) then Ci + C2 + Cs + C4 = and 
for each i E {1, 2, 3, 4} we have 



Ai = {xeR-': {x, Ci) = max {x, (j) } ■ (6) 

t i6{l,2,3,4} J 

By [23l Lem. 3.3 & Cor. 3.4] the vectors {Ci}i=i are distinct, nonzero and not coplanar. 

Since the sets in question are cones, it is beneficial to study them in terms of their intersection 
with the sphere S'^, i.e., define 

Ti = AinS^, (7) 
and letting a denote the surface area measure on 52 (thus cj(52) = 47r), define 



xda{x) = V27r3 • (8) 

IT, 

where the last equality in ([8]) follows from integration in polar coordinates. Being a constant 
multiple of the vectors {Ci}i=i, the vectors {zi}f^i are also distinct, nonzero, not coplanar, and 
satisfy zi + Z2 + + Z/^ = 0. Moreover, 

Ti i X £ : (x,Zi) = max {x,Zj) \ . (9) 

[ i6{l,2,3,4} J 

Consequently, 

Vi,j e {1,2,3,4}, xeTiHTj ^ {zi,x) = {zj,x). (10) 

Fix £ € {1, 2, 3, 4}. Since {zi}f^^ are not coplanar there exists a unique V£ G 5^ satisfying 

V«,j E {1,2,3,4} \ {^}, {zi,ve) = {zj,ve) > {ze,ve). (11) 

The vectors {ve}j^i are the vertices of the partition P = {Ti, T2, T3, T4} of S*^. A simple argument 
presented in Section [3] shows that Vi ^ {—Vj,Vj} if i ^ j and each spherical triangle Tj is contained 
in an open hemisphere of 5^. Moreover, for distinct G {1,2,3,4} we have det{vi,Vj,ve) 7^ 0, 
and if we define 

9ii = arccos((fi, f ,)) and Qijg = arccos ( ( -rr— 77-, rr^ 17- )) > (12) 

\\\\vi X Vj\\2 \\Vj X Vl\\2/ J 
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then it is also argued in Section [3] that 9ij,@iji S (0,7r). Observe 9ij = 9ji, @ijk = @kji, and Qiji 
is the spherical angle, at vertex Vj, of the spherical triangle with vertices {vi,Vj,Vi}; the cosine of 
this angle is exactly the inner product of unit normals of the two planes containing {vi,Vj,0} and 
{vj,Vi,0} respectively. See Figure [2] for a schematic description of the notation. 




(a) (b) 



Figure 2. (a) Spherical triangle constructed with three great circles (b) Partition 
of the sphere into four spherical triangles. 



In order to proceed we need to have a formula for the spherical moment Z4 in terms of the vertices 
{vi,V2,V3} of the spherical triangle T4. In order to do so, assume that det(i)i, 1)2, t's) > 0. Then 

1 VtXV?. V2XV3 V3XV1 \ 

+ '^2377— ——71 — r "3177— ——rr ) • vJ-"jj 



Z4 



'121 



\Vl X V2\\2 \\V2 X Vs\\2 ' I I'^^S X "^1 1 12 / 

Identity ([T3D was proved by Minchin in 1877; see [321 p. 259]. We also found identity in [S], 
which is a military publication that is not publicly available (it isn't classified: we purchased 
access to it). Since both references for are hard to find, we include a brief derivation of it in 
Proposition 13.21 using the Cauchy projection formula (see j26|. p. 25]). 

Using (jl3|) and geometric arguments, we show in Lemma [3.51 that the vertex is determined by 
the vertices {vi,V2,V3} as follows; an analogous formula expresses each vertex of the partition in 
terms of the other three vertices. 

1 

V4 



{ sin 023 sin0i3 sin0i2 
' ""1 H 7, V2 H -p. 



VX V ^23 



where 
A = 



sm ^23 sm t/i3 sm U12 
—p, vi H V2 H U3 



723 



713 f7i2 

sin 6*23 sin 6*13 



713 



sin^ ^23 sin^ 

I "Tr 



712 



713 



02 
^723 



sm fc'12 
'712 



^, ,, „ smt/23Sinf7i2 „ smt/i2smt/i3 „ , , 

+ 2 f^^T COS 012 + 2 fi- COS 013 + 2 COS 023. (14) 

f723fc'l3 fc'23'712 t/i2t7l3 

In Proposition 13.71 we derive the following restriction on products of opposite edges. 

sin 012 sin 034 sin 0i3 sin 024 sin 0i4 sin 023 



7l2t734 



7131724 



7141723 
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Also, in Lemma 13.91 we obtain the following expression for the objective function at the partition 
P = {Ti}f^^ as the sum of the squares of the angles between the vertices. 



Define 



FiP) = Ikilli + \\z2f2 + IWsWl + Ik4||i = + oh + + Bl^ + Ql^ + Ql^. (15) 

sin 6^23 . , sin 6*13 sin 6^12 „ 

7 = — B cos 6*12 + —T \ cos6'23. (16) 

C'23 f7i3 6^12 

In Lemma 13.81 we obtain the following useful restriction on a single spherical triangle appearing in 
the partition {Ti, T2, T-^^T/^. 



cos 



^^ sin£!Lj = ~' (V»;. + «i3 + 2«,,fecose,3) = (17) 



where A is given in ([14 

Observe that (jl7p is a restriction involving data from only one triangle T4, as we see by the 
definitions of 7 and A. Moreover, by symmetry, every cyclic permutation of {^125^237^13} must 
also satisfy (jl7p . And, every triangle of P must satisfy (|17p . with the indices of the Qij substituted 
appropriately. To understand the ramifications of (fTTl) assume that 6*12 = ^23 = ^13 ='■ d- Then 
(fTTl) simplifies to 

cos f \/2(2cos6' + l)(l -cos6l)-^ ) = — ^\/2cos6l + l. (18) 
V sm fc* y V3 

The only (9 G (0,27r/3) satisfying this equation are 6 = arccos(-l/3) ^ 1.9106332362490187, and 
9 ~ 1.5379684120790425. The first solution corresponds to the partition whose vertices are those 
of a regular simplex inscribed in the sphere, and the second intuitively corresponds to a critical 
point "between" the regular partition and the propeller partition (note that 9 = 2tt/3 implies 
det(fi, f2, ^3) = 0, and 9 > 27r/3 cannot be achieved by a spherical triangle). The above special 
values of 9 will play an important role in the ensuing computations. 

As explained in the discussion following Lemma [3.91 ^ combination of (I15p and ()17p can be used 
to express the value of our objective function at the maximizing partition P in terms of the data 
from a single spherical triangle in P. Specifically, we have 

F{P) = 3 (0^2 + Oh + Oh) + 2 COs(e213)^12^13 + 2 COs(ei23)ei2e23 + 2 COs(e23l)^23^13. (19) 



In particular, it follows that the right hand side of ()19p must be the same for each set of data for 
all four triangles of P. 

Up to now we only discussed identities that the extremal partition P must satisfy. In order to 
proceed we need to prove some a priori estimates on the various parameters in question. First, let 
M = max{6'i2, 6*13, 6*23}. It follows from ^ that F{P) ^ UM'^ (using the fact that at least one 
@ijk must be greater than ^). Our contrapositive assumption ([5]) combined with ([8]) implies that 
F{P) > so we deduce that all the triangles in P must have an edge of length greater than 
^'^ > J- Additional arguments (that are significantly more involved technically) imply that for all 



2V14 

distinct i,j,£ G {1,2,3,4} we have 9ij ^ vr — ^ (Lemma 13. 12p and 9ij,smQij£ > ^ (Lemma 13. 16p . 

In Lemma [4.11 we prove that (^12, ^13, ^23) must be outside the £2 ball of radius ^ centered at 
{9, 9, 9) for 9 = arccos (-i) and 9 = 1.53796841207904. This excludes from our search space balls 
centered at the two solutions of (jl8p that were discussed above. In Lemma 14.21 we show that, up 
to a relabeling of the spherical triangles Ti,T2,T^,T4^, we may assume that \/A > 

It turns out that the above estimates suffice in order to conclude the proof of Theorem 11.11 via 
a search over a sufficiently fine net. To explain this endgame, note that by the spherical law of 
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cosines we have cos ^13 = cos Ou cos 6*23 + sin 6*12 sin 6*23 cos 6123 • It therefore follows from plj) . ([T 
(fTHj) that if we define /i, A, 7 : [0, vr]^ ^ M by 

sin^ X sin^ y sin^ z sin x sin y sin x sin z sin y sin z 

A(x, y, zj = 1 1 h 2 cos z + 2 cos y + 2 cos x, 

a; y 2 xy xz yz 

sinz siny sinx 
7(x, y, zj = cos X H 1 cos z, 



7 I \ m r / / o o ^ cos y — cos X cos z \ , , 

h[x,y,z) = VA(x,y,z)cos Wx^ + z^ + 2x2; +7(x,y,2;), 



sm X sm z 



then h{6i2, ^13, ^23) = 0. This identity must hold for all cyclic permutations of the indices {1, 2, 3}, 
so we also have the two identities /i(^i3, ^23i ^12) = ^(^28)^12)^13) = 0. It follows that if we 
define H = {Hi, H2, H3) : [0, vr]^ — )• by H{x, y, z) = {h{y, z, x), h{x, y, z), h{z, x, y)) then at the 
extremal partition P we have H{6i2, Oi^, 623) = 0. 

Our strategy is therefore as follows. At every point q in the search space given by the constraints 
described above (specifically, see the system of equations (f67|) . excluding also the two balls described 
in Lemma l4.ip . we will show that either H[q) 7^ or Fq^q) < 97r^/4, where Fq is given by the right 
hand side of (fT9|) (with the understanding that one expresses cos 6123 , cos 6213 ) cos 6231 in terms 
of 9i2, ^23) 6*13 using the spherical law of cosines, so that Fq is a function of 6*12) ^23) ^13)- Due to the 
above discussion, such an assertion will complete the proof of Theorem ll.il Moreover, it turns out 
that this assertion holds "with room to spare" , making a computer-assisted verification feasible. 

To this end, we need to have estimates on the modulus of continuity of H and Fq. Such esti- 
mates are complicated but can be proved using elementary considerations; see Lemma 14.31 and the 
discussion immediately following it. Consequently, if we are given a point q in our restricted search 
space and r G (0, 1) for which either \H{q)\ > r or FQ{q) < 97r^/4 — r, then it would follow that for 
some r > that depends on r via our modulus of continuity estimates, in the entire ball of radius 
r centered at q either H ^ or Fq < 97r^/4. 

From this procedure we achieve a process by which we can iteratively remove from our search 
space macroscopically large regions of controlled size in which we are guaranteed that a coun- 
terexample to the Propeller Conjecture cannot exist. Our program proceeds to remove such balls 
until it eventually exhausts the entire search space, thus arriving at the conclusion that there is no 
counterexample to the Propeller Conjecture. 

In order to make such a procedure rigorous, we carefully account for all possible rounding errors. 
We do this quite conservatively, i.e., while significantly overestimating the magnitude of all possible 
errors, as explained in detail in Section Note that our computer-assisted proof does not compute 
any computationally complicated expressions such as, say, definite integrals or implicitly defined 
functions: it only checks inequalities between expressions involving elementary combinations of 
trigonometric functions and square roots. 

In summary, the validity of the Propeller Conjecture in has mainly conceptual implications, 
and what remains open seems to be of a more technical nature: to perhaps find a clever transfor- 
mation, e.g., as in Ball's cube slicing theorem [3], that allows one to prove the theorem analytically 
rather than resorting to a transversal of a net in a region where the desired inequality is actually 
quite weak. At present, we have a clean geometric argument which addresses directly the region 
where the Propeller Conjecture is most subtle, and in the remainder of the search space we do 
something rather crude. This has two conceptual consequences: the Propeller Conjecture as a 
nonintuitive link between complexity theory and geometry is correct, and moreover one has a geo- 
metric/analytic proof of this conjecture in a region where it is tightest. Nevertheless, we are lacking 
a technical idea that allows us to address the remaining region in a way that does not resort to 
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"brute force" . It remains a challenge to find such an idea, with the hope that it will pave the way 
to a proof of the Propeller Conjecture in all dimensions. 



This section contains proofs the geometric identities and inequalities that were stated in Section^ 
Recall that we are assuming that P = {Ti, T2, T3, is a partition of S*^ into four spherical 
triangles, and that {zi}^^^, as defined in ([5D, are the corresponding spherical moments. We are also 
making the contrapositive assumption that our objective function 



(recall (jlSp ) exceeds 97r^/4, which is the maximal value of F{-) that is predicted by the Propeller 
Conjecture. We also assume that P maximizes F. The vertices of the partition {vi}f^i were defined 
in Section [21 along with the angles Oij and Qiji as given in ()12p . 

Proposition 3.1. (a) zi + Z2 + + = 0. 

(b) Suppose two triangles Ti,Tj € P share an edge E = TiDTj. Then for all e £ E we have 
{zi,e) = {zj,e). In particular, if i ^ {hj} then {zi,ve) = {zj,ve). 

(c) For all distinct i,j E {1,2,3,4} we have {zi,Vi) = {—3zj,Vi). 

(d) For all distinct i,j S {1,2,3,4} we have Vj ^ {—Vi,Vi}. 

(e) For all distinct i,j,i G {1,2,3,4} we have det{vi,Vj,ve) / , each Ti is contained in an 
open hemisphere, and < Oij, Qijg < vr. 

Proof. Parts (a) and (b) were already proved in Section [2l for part (b) see (jlOp . A combination 
of (jlip and part (a) implies part (c). For i £ {1,2,3,4} let II^ be the unique plane containing 
{-^i}'t6{i,2,3,4}\{^}- A vector is perpendicular to II^ if and only if it has equal inner product on 
{-2^i}'t6{i,2,3,4}\{^}) so Wv£ is also perpendicular to Hi. Since {zi}f^^ are not coplanar the planes 
{n^}^^-^ have distinct perpendiculars. So vi,V2,vs,V4 are distinct with Vi ^ M.Vj for i / j. This 
establishes (d). For every i G {1,2,3,4} consider the set 



It follows from (jlip that Ug is an open neighborhood of vi, and in combination with ([9|) we know 
that 



ie{l,2,3,4} \{£} ^ UinTi = IxeUgDS'^ : {x,Zi) = max {x,Zj)\. (21) 



Write / = {1, 2, 3, 4} \ {i}. For every distinct i,k £ {1, 2, 3, 4} the set {x £ S'^ : {x, Zi) ^ (x, Zj)} is 
a half sphere, and therefore the set Ti^i = {x £ S'^ : [x, Zi) = maxjg/(x, Zj)} is strictly contained in 
a half sphere (since the Zj are distinct and nonzero). Note that dTi^i is the union of two half circles 
that meet at antipodal points, and vi is one of these antipodal points. So, the spherical angle of 
Ti^i at the vertex vi is in (0, vr). By ([9]) and (j2ip this spherical angle is identical for Ti^i and Tj. This 
shows that G (0,7r) for all distinct i,j,l G {1,2,3,4}. Also 6ij G (0, vr) since Vj ^ {—Vi,Vi}. 
Next, assume for the sake of contradiction that <lei{vi,Vj,V() = 0. Since Vi ^ Mfj for distinct 
i,j G {1, 2, 3, 4}, it follows that Vi = ajVj + aivi with aj, G M \ {0}. Hence Vi x Vj = a^v^ x Vj 
with ai 7^ 0, contradicting Qiji G (0,-7r). We have already seen that Ti^i is contained a closed 
hemisphere, and this containment can be chosen such that Tj^^ intersected with the boundary of 
the hemisphere is {v^, —V(}. But from @ we know that T, C Tj ^ and Tj avoids a neighborhood of 
{—Vi}, so each Tj is contained in an open hemisphere. This concludes the proof of (e). □ 



3. Proofs of basic identities and estimates 



F{p) = Il-2illi + lk2i + lk3lli + 1^4111 



(20) 
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We now give an explicit formula for the spherical moments {zj}^^]^; as explained in Section [21 
this formula can be found in j32|, p. 259] and [9l p. 6]. Since these references are hard to find, we 
provide a proof. 

Proposition 3.2. Suppose that the spherical triangle C S"^ has vertices {vi,V2,vz^ satisfying 
det(ui, U2, U3) > 0. Then the vector Z4 = xda € satisfies 

U = \{ ^12/^ + ^23/^ + ^31/^) . (22) 

2V II'^1X^2||2 \\V2XVz\\^ IhsXUillg/ 



Consequently 



I V / -^23 ^"2 X \ 6l23det(z;i,W2,W3) ,oo^ 

{zi,vi) = {—- —^vi) = — — (23) 

^ 2 i;2 X ■i;3 2 / 2sm6'23 



Proof. Fix y G r4 \ dTi^. By definition of Z4, 

(^4,2/) = / {x,y)dc7{x). (24) 

Let C ]-,g i^jjg convex hull of U {0} and let {C/i}f=i ^ be half spaces through the origin 
such that {x G : ||x||2 1} fl (flLi = H; see Figured 




Figure 3. Spherical triangle together with three disc sectors, each containing {vi,Vj,0}. 

Define H: 5^ — )■ so that n(x) is the orthogonal projection of x S S*^ C onto the unique 
plane that is intersecting the origin and perpendicular to y. Applying the Cauchy projection formula 
(see for example [261 P- 25]) to T4, we see that 

/ {x,y)da{x) = AreaiR2n(r4 n {x G S"^ : {y,x) > 0}) - AreaR2n(r4 n {x G S'^ : {y,x) ^ 0}). (25) 

Define A = {T4 \ 8X4) n {x G : {y, x) > 0} and B = [OH) \ A. Note that dH is the union of 
T4 and three sectors of discs {Di}^^^ such that Di C dUi. Since y G C\i=i the exterior normal 
n(x) of X & Di CI dH has nonpositive inner product with y for i = 1,2,3. So, {n{x),y) > for 
X G dH if and only if x G A. Since H is convex, H: A ^ Il{dH) and IT: S — t- Il{dH) are bijections 
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almost everywhere. Since dH = T4 |J ^U?=i ^ij > the Cauchy projection formula applied to right 
side of the identity Aiea.^2{Jl{A)) = (J1{B)) gives 

AreaK2n(r4 n {x G 5^ : (y, x) > 0}) = Area]R2n(r4 n {x G 5^ : {y, x) ^ 0}) 

6*12 / Wl X 772 \ , 6^23 / ^^2 X \ 6*13 / ^3 X t;i \ (26) 

+ ^ ( II II ) + ^ ( II II ) + ^ ( y, ■ 



2 \ X 772112 / 2 \ ||i;2 X t;3||2 / 2 \ 11^3 X uill 

Since ^ holds for any y G n r4 , ([MD , (ES]) and ([26]) give ([22]). □ 

Note that the vectors in the determinant of ()23p are assumed to have positive orientation. So, 
when we apply this equation below we will need to respect orientations. 

Corollary 3.3. Under the assumptions of Proposition [37B. we have 

Ik4||i = ^(e?2 + ^i3+^?3-2cOs(e213)^12^13-2cOs(ei23)^12^23-2cOs(G23l)^23M (27) 

- ^ fff2 ifl2 ,n2 , ^ COsgi2COsgi3-COsg23 ^ . 

- 4 ^^12 + ^23 + ^31 + 2 —ff—^^ ^12^13 

COS 012 COS 6^23 - COS 6*13 ,„ COS 023 COS COS 6*12 „ „ \ 

+^ ^-7; ^-7 yi2c'23 + ^ ^-p. cy23c'l3 • V^oj 

sm O12 sm 023 sm 023 sm 0i3 J 

Proof. Apply (j22p . and then use either the vector identity {a x b,c x d) = (a, c){b, d) — {a, d){b, d), 
or the spherical law of cosines. □ 

Corollary 3.4. For distinct i,j G {1,2,3,4} we have {zi,Vj) > 0. 

Proof. Let k,i be such that det{vk,Vj,V£) > 0, and apply (j23p and Proposition 13. iT e). □ 

Lemma 3.5. Write {z4,vi) =: a, (2:4, ^2) =• ^> (-24, ^"3) =: c and (24, ^4) =: —3d. Then the following 
relation holds 

1111 

-Vi + -U2 + -V3 + -V4 = (29) 

a c a 

Proof. By Proposition I3.1f e) we have det(fi, V2,V3) ^ 0, so we may write = X]i=i ^i^i some 
ai,a2,a3 G M. By Proposition 13. iT c). 

d= {zi,V4) = (2:2, 774) = (23,i;4). 

Substituting in the expression = Yl^=i ^i'^i then applying Proposition 13.1( b) and (c) several 
times, we get 

d = — 3aai + 602 + ca3 = aai — 36a2 + ca3 = aai + ha2 — 3ca3. 
Solving this system of equations gives 

d = —aai = —ha2 = —ca^. 
By Corollarv 13.41 we may write ai = —d/a,a2 = —d/h,a^ = —d/c. Substituting these equalities 
into the expression 774 = Y^=i ^i'^i gives (|29p . Note that at least one ai is nonzero, so d is nonzero 
as well, that is, we have not made a division by zero. □ 

Since det(t;i, 772, ^3) / 0, we have \\^vi + j;V2 + ^^3112 / 0. Since 11774112 = 1, Lemma [331 implies 

^4 - -771 -T -T — rp- l-^Uj 

fe?^2 + ^^^3||2 

A priori, we have only determined 774 up to multiplication by ±1. However, Proposition 13. iT c) and 
Corollarv 13.41 implv that d = (—1/3) (24, 774) > 0. So, taking the inner product of Z4 with both sides 
of ipU]) shows that ([50]) has the correct sign. 
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Applying ([23]) to ([SOD gives 

f 4 = 

Define, as we did in Section [21 

sin 6*23 



sine23 
023 



ei3 



sin 023 
623 



-1 + ^-2 + ^ 



V3 



A 



smt/is sint/i2 
^1 H 7 "^2 H 7, ^"3 



723 



713 



712 



sin2 6'23 , sm^0i3 sin^ 6*12 , sin ^23 sin 6*13 

H — \ h 2 — cos 012 



92 

^23 



^713 



'712 



723f7l3 



sm 6*23 sm tyi2 „ , sm t/12 sm ^13 
+ 2 — cos 1^13 + 2 — cos 6723. 



(31) 

(32) 
(33) 



^23^12 ^12^13 

Proposition 3.6. (a) All four vertices of P are not contained in any closed half sphere. 
(b) The function F{P) is well-defined as a function of the vertices {vi}f^i. 

Proof, (a) By Proposition 13. li fe) we have < 6ij < ir. (j3ip says that ^4 is in the complement of all 
closed half spheres containing {vi,V2,vs}. 

(b) From Proposition 13. li fe), the edges Ti n Tj of P are geodesies of length less than vr, so the 
vertices uniquely determine the edges. Finally, the edges determine P itself. □ 

Proposition 3.7 (Restrictions on products of opposite edges). 

sin 012 sin 6*34 _ sin ^13 sin 024 _ sin ^14 sin ^23 ^^^-^ 

^12^34 013024 0U023 

Proof. From Proposition 13.1( c) we have {z3,vi) = {z2,vi). For each side of the equality {z3,vi) = 
{z2,vi), apply ([23]) for {zi,vi) and then substitute in (f3T]) to get 

det(t;i,t;3,7;2) sin 6*12 6*24 _ det(fi, t^a, ^2) sin6>i3 6*34 
2y/X 012 sin 6*24 2\/A O13 sin 6134' 

Canceling terms gives the first equality of (j34p . The second follows similarly: Proposition I3.1f c) 

says {z3,V2) = {zi,V2), and so on. □ 

Lemma 3.8 (Restrictions on a single triangle). Let A be defined as in (j33p and let 7 be defined as 
in ([TB]) . i.e., 



smy23 o , smt/i3 sin«i2 „ 
7 = —p. cos 012 H p 1 p cos 6^23- 



Then the following holds 



723 



cos 



713 



7231712 



712 



sm ^23 sm 6'i2 



7 



COS 



2^ + 02^ + 2012^23 COS 6123 ) = 



(35) 

(36) 
(37) 



Proof. Arguing as in Proposition 13.71 we start with Proposition I3.1f c) and the equality (24,^1 
(23,^1). For each side of this equality, apply (j23|) for {zi,vi) and substitute in (jSTj) to get 

^23^12 ^24 1 



sin ^23 sin 0i2 sin 6*24 \/A 
Substituting ^24 = arccos((i;2> '^'4)) into ([38]) and rearranging gives 

^23^12 



(38) 



cos 



^\/Aa/1 - {v2,V4y' 



sm ^23 sm t/i2 



{V2,V4). 
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From ([3T|) and the definitions of A and 7 ( (p3]) and (j35]) ). we have {v2,V4) = —j/^/X, so 



cos 



7^ 6'236'l2 \ 7 



A sin 6*23 sin 6*12 J ^/x' 
yielding ([36j) . To derive ([37|) . observe that A — 7^ allows some cancelation as follows. 



.2, 



2 sm ^23 2 /I \ sm t/12 2 /I N sm t/23 sm t/12 
A-7 = — ^2 — (1-cos di2)-\ T2 — (1 - cos 6/23) + 2 ^— T (cos6'3i - COs6'23COs6'i2) 

^23 ^12 c/23t/i2 



sin^ 6*23 ^ • 2z) A , 012 , . 2/) N , ^ sin6'23 sin 012 ^ /j /j n 

= — To — (sm 6^12) H T2 — (sm 6^23) + 2 r (cos t'si - cos 6^23 cos 6^12). 

^23 (712 t723t7l2 

Thus, 

(\ ^,2N^_^|3^1^_ fl2 , zj2 , 9/1 cos 6*31 - cos gas cos gi2 
J— . 2/1 = ^12 + ^23 + 2fc'l2t'23 ^"7 

sm^ ^23 sm^ Q\2 sm 6'i2 sm 6*23 

= ^12 + ^23 + 2^12^23 COS ©123 , by the spherical law of cosines. 



Substituting this equation into (1361) gives (j37l) . □ 
Since — = {v4„V2) = cosg24, (|37|) says that 

^?2 + ^i3 + 2^12^23 COS Gi23 = 0^^. (39) 

Substituting (j39p back into (|27p and (j20p gives a simplified form of F{P). 
Lemma 3.9. 

F{P) =5^4 = ^(arccos((^,,^,)))2 (40) 



^4 II ||2 



Proof. Applying (j28|) to the definition -F(P) = Yli=i \ \^i\\2 gives 

^(^) = ^ E 4 - ^ E cos(e,,-,)%g,-,. (41) 

Here the first sum in (j4ip runs over all (i,j) satisfying 1 ^ i < j ^ 4, and the second sum runs 
over all equivalence classes of 3-tuples under the equivalence ~ where 

are distinct elements of {1,2,3,4}. Prom ([39]) we have, for r ^ 

el + + lOijO,^ cos e,jf = (42) 

Substituting into dH]) gives 



^(^) = ^E^?,-7E(^J^-4-^l.)- 



2^ u 4 

In the right summation, for a fixed with i' 7^ j' , the term g^,^ appears four times with a 

minus sign and two times with a plus sign. Therefore, 

^(^) = ^E^?.-iE(^l^-4-^l.) = ^E4-i(E-2^?,) = E4' 

proving (fiOl) . □ 
In (HO]) replace 6l|^ with the left side of ([12]) for (j,r) G {(1,4), (2,4), (3,4)} to get 

F{P) = 3 (0?2 + ^23 + ^Iz) + 2cOs(e213)gl2% + 2 COs(ei23)gl2g23 + 2cOs(e231 )6'236'l3- (43) 
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Definition 3.10. By the spherical law of cosines, Qijg is a function of {9ij,9ji,6ii}. Let 

Fo{ei2, 023, On) = 3 (0?2 + Oh + Oh) + 2 COs(e213)^12^13 + 2 COs(ei23)^12^23 + 2 COs(e231 )023^13. 

By (j43p . F{P) = Fq can be computed by the data from only one triangle of P. That is, 

-^0(^12,^^23)^13) = -^0(^12,^24,^14) = -^0(^34,^24,^23) = -^0(^34,^14,^13)- 

Since the spherical propeller partition, i.e., the partition of S'^ corresponding to the intersection 
of the propeller partition of M"^ with S'^, satisfies F = (9/4)7r^, we may bound (I43p in terms 
of M = max{^i2, ^13, ^23}- In particular, we get F{P) ^ 14Af^, which is less than (9/4)7r^ if 
M < TT (In our bound on F{P), we use that one @ijk must be larger than 7r/3 by the 

pigeonhole principle applied to 6123 + ©231 + ©312 — vr = Area52(T4) > 0, so one cosine term 
satisfies cos0jjfc ^ 1/2.) So, this bound on F gives 

Corollary 3.11. Any triangle Ti of P must have an edge of length greater than ir > |. 



'^j<^^-n■ (44) 



Iteratively applying ()34p gives the following improvement to Proposition IS.lT e) . This improve- 
ment will be important in our numerical calculations in Section [4] and in the proof of Theorem 11.1 
since we eventually require a bound on the derivative of 9ij/ sin 9ij. 

Lemma 3.12 (Restrictions on edge length). For every distinct i,j £ {1,2,3,4}, 

1 
2 

Proof. Suppose 9i2 ^ vr — 1/2. We are eventually going to derive a contradiction, which will 
allow us to conclude that 9i2 < vr — 1/2. We will essentially only need (j34p . the monotonicity of 
X (sina;)/x for x £ (0, vr), the fact that the perimeter 9ij + 9ji + 9£i of a spherical triangle is 
bounded by 2vr, and that (sinx)/x ^ 1 on [0, vr]. Relabeling the edges if necessary, we may assume 
that 6*13 = max{6'i3, 6*24} • Using 6*12 ^ vr - 1/2, dS]) gives 

sin(vr - 1/2) ^ sin 6112 ^ sin 6*12 sin 6134 _ sin 6113 sin 6*24 ^ / sin^y ^^^^ 



(46) 



vr — 1/2 9i2 ^12^34 C'13'724 V Cl3 

Therefore, 6*13 ^ X where X E (0, vr) satisfies 

sinX\^ sin(vr — 1/2) 



X J vr-1/2 

So, 013 ^ X > 2.065. Now, since 9i2 + ^13 + ^23 ^ 2vr the bounds on 9i2 and ^13 give 

923 ^ 2^ - (vr - 1/2 + X)=:Y < 1.577. 
But then (p4l) and the bounds on 9i2 and ^23 yield 

sin0i4 sin0i2sin034 ^23 ^ sin(vr— 1/2) Y 



Oi4 ^126*34 sin 023 vr — 1/2 sinY' 

So, 014 ^ Z where Z G (0, vr) is defined by 

sinZ _ sin(vr - 1/2) Y 
Z ~ vr-1/2 siny' 
Then ^14 ^ Z > 2.388. But since 0i4 + 024 + ^12 ^ 2vr the bounds on 0i4 and 0i2 give 

024 ^ 2vr - (vr - 1/2 + Z) =: W < 1.254 
We now iterate the above argument. Using (|45p . 



sin 013 _ sin 012 sin 034 024 ^ sin(vr - 1/2) W ^^^^ 



?13 C12C34 



sin 024 vr — 1/2 sinp^ 
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So, 6*13 ^ X with (sinX)/X = (sin(7r - l/2)/(7r - 1/2))(VF/ sin T^). We now repeat the argument 
above, starting at (fl6|) and using tildes to designate updated variables. Thus, ^13 ^ X > 2.499, 
6*23 ^ < 1.14,6*14 ^ Z > 2.52,6*24 ^ # < 1.122. Using W in place of W in gT]), we find that 
% > 2.52. 

In summary, iterating the above procedure improves our bounds on the Oij. We find that 
6'i3,6'i4 > 2.52. We also find that 6*23, 6*24 < 1-14. Since 6*13 + 6*14 + 6*34 ^ 27r, 6*34 ^ 27r - (2(2.52)). 
But then 6*23, 6*24, 6*34 < 5/4 < 71—4=, violating Corollarv 13. 1 II and giving our desired contradiction. 



2V14' 

We therefore conclude that Oij < vr — 1/2. □ 

For a partition P = \Ti\\^Y ^ above, the following map G: (5"^)^ — )■ M is also a function of P 
by Proposition I3.6l fb) 

G'(i;i,^;2,^;3,^^4)= 4- (^8) 

Proposition 3.13. Suppose that the conclusion of Proposition [WHY e) holds. Then Proposition 
\3.lV c) holds if and only if P is a critical point of G. 



Proof. For a manifold M with v G M, T^M denotes the tangent space of M at v. Since G = 
^i<j^ij ~ X]i<j(^^*^cos(fi, I'j))^, we take derivatives of G with respect to Vi. For ^2,^3,^4 G S"^ 
fixed and x G variable. 



VS2G{-,V2,V3,V4) = ( ViRsG ( ■n-^,U2,U3,U4 



5-2 



Recall that T^j^M^ is isomorphic to M'^. We denote this isomorphism by [d/dy] -H- y. 
Let ^ be a vector in the tangent space T^^M^. Identifying ^ with y G M"^ we con 

d^f X \ v-^ — 2arccos((x,tio)/ llxIL) / d, .,, ,, , 

dy Vlkll2 / ^ y'l-{{x,Vj)/\\x\\^y \dy 

^ -2arccos((x, Vj)/ ||x||2) f 11^112 ^(^'^i) ~ ^(^'^i) 11^112 ^ ^(^'^) 



\/l - ((2;, ^^i)/ 1^112)^ V ll^llz 

— 2arccos((x,t'j)/ IIXII2) / / ||x||2 y — x ||x||2 (y, x) \ 
^l-{{x,v,)/\\x\\,r \ I 1 ' 



Letting ^ G T^^S"^ and x = f 1 G 5^ gives (y,x) = 0, so that 



d / X \ ^ -2arccos((?;i,7;j)) -r-^ -26'ij 

7^-r^,t'2,t'3,'f4 = > ^======^ = > ^— — y,?; ) . 

So, if we take ^ = x t)2 G T^^S"^ and then use (p3|) . we get 

■^G = .^^^^ (lei{vi,V2,V'i) ^^^det(-yi,f2,';^4) = 2(-(z4,^;2) + {z'i,V2)). 

dy sm ^13 sm 6^ 

In general, if we take the derivative at vi in the direction v-i x Vj and set this derivative to zero, 
we get {zr,Vj) = {z£,Vj) where the set {i,j,£,r} is equal to the set {1,2,3,4}. Since the above 
argument can be reversed, the proposition is proven. □ 

Setting ^ = f 2 — {v2,vi)vi and -jj^ = V3 — {v3,vi)v3 in the above argument gives 

- 6*12 = 6*13 cos 0213 + ^'i4 cos 6214- (49) 
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Observe that 

( Trir^'"2,V3,vA = V -T^^ {y,Vj) = V . ^^^^ ((^2 - {vi,V2)vi),Vj) 
= - ((1^2, t'l))^) - -?^^{{V2,V3) - {{v2,Vi)){{vi,V3))) 

sm (712 sm (713 

- ■ {{V2,V4) - {{v2,Vi)){{vi,V4))) 

sm (714 



-(1 - {{V2, vi)) ) - . . ((til X t;2), (vi X V3)) - . „ ((ui X z;2), (I'l x i;4)) 



sm (^12 sm Wi3 sm ^14 

Thus, dividing by sin0i2, setting = 0, and simphfying as in Corollary 13.31 gives (j49]) . 

We now give a geometric interpretation of (1390 and ()49p that should clarify their meaning. Define 
expfi; {2,3,4} by 

expr'(i)=0.. ,^''°""f;' . (50) 

Ibi - W''^i)^'ill2 

Then exp^^(j) is perpendicular to "Wi and parallel to Vj — Proj^^fj, with length 9ij. 

Lemma 3.14. Suppose that the conclusion of Proposition [SHV e) holds. If {vi}f^i is a critical point 
of G, then 

exp^i(2) + exp^i(3) + exp^i(4) = 0. (51) 

Proof. Let Yiij C be the plane containing {0,i)j,fj}. For j G {2,3,4}, exp|f^(j) € Ilij is 
perpendicular to vi and in Hij, so (exp]~^(2), exp^^(3)) = ^12^13 cos 0213- We can therefore rewrite 
(f39D as 

||expr'(2)+expr'(3)||2= ||exp^i(4)||2- (52) 
Using the definition of exp]"^(j), we similarly rewrite (j49p (with its indices permuted) as 

- expr^(4) = Proj^^p-i(,)(exprH2)) + Proj^^p-i(,)(expr^(3)). (53) 

Since exp^^(j) is perpendicular to vi for j G {2,3,4}, {exp|f ^(j)}jg|2^3^4} lies in a plane containing 
the origin. Combining (j52p and (j53p therefore gives (|5ip . □ 



In particular, orthogonally projecting (j5ip orthogonal to exp^^ ^(2) gives 

6*13 sin 6213 = 6*14 sin 6214- (54) 

The similarity of ()39p . (j49p and ()54p is no coincidence. For example, note that (j49p and (j54p 
express X]jg{2 3 4} 6xp^^(j) = 0, tangential and perpendicular to exp]~^(2), respectively. So, com- 
bining Proposition 13.13] (j49p . (j54p and Lemma 13.141 we deduce the equivalence 

Corollary 3.15. Suppose that the conclusion of Proposition [0] (^e ) holds. Then 

{Prop. \3l\fc)\ ..Pisa i fdSH) /^o/fis at eac/i Vi. T/iai is, 

/ {critical point of Gj ]^for each i, -^^ exp. ^(j) = 

Lemma 3.16. O-ij > 1/10 andshiQij^ > 1/10. 

Proof. After recording a few identities, we consider a few cases in proving Oij > 1/10. In all cases, 
we argue by contradiction and assume < ^?i2 ^ 1/10. Also, by relabeling if necessary, we may 
assume that max{0i3, ^23) ^14> ^24} = max{0i3, ^23}- 

cos 6*23 - cos 6'i2 cos 6*13 

cos 0213 = r-7 ■ (55) 

smcyi2 sincyi3 
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COS 6*12 - COS 6*13 COS 6*23 

COS 0132 = — — . (56 

sm 6^13 sm 6^23 

COS 6*34 = COS ^13 COS ^14 + sill ^13 sill 6*14 COS G314 

(57) 

= COS 6*13 COS 014 + sin 6*13 sin 6*14 cos(27r — G213 — 0214)- 
I arccos(x + ?/) — arccos(x)| ^ arccos(l — y), x,y,x + y £ [—1,1]. (58) 
(Case 1) 7r/2 ^ max{0i3, ^23} ^ 7r/2 + 1/10. Without loss of generality, label the edges so that 
^13 ^ ^23) so 7r/2 ^ ^13 ^ 7r/2 + 1/10. Then 614 ^ 624, from Proposition I3.1l fb). The Proposition 
shows that Projgpan{?;i,t;2}(-^4~-2^3) = 0- Suppose On < 624- Let e be the midpoint of the edge r4nT3 
and let n G span{ui,z;2} be perpendicular to e such that {n,V2) > 0. Since ^13 ^ 623, {z4,n) ^ 0. 
To see this, let 11: 5^ — t- denote projection onto the unique plane intersecting the origin and 
perpendicular to n. Without loss of generality, the vertices are oriented so that det('yi, 112) ''^3) > 0. 
Let m be the point of intersection of this plane with {dT^) \ {dT^). The spherical triangle TJ^ 
with vertices {vi,V2,m} then has two equal length edges {vi,m} and {v2,m}, and T4 C T4. Also, 
6^13 ^ 6*23 implies that T^nixG S^: {x, n) ^ 0} = n {x G S'^ : {x, n) 0}. Since T4 D and 
is isosceles, 

n (T4 n{xes^: {x, n) ^ o}) ^ n (r^ n {x g 5^ : (x, n) ^ o}) = n [T^ n {x G 5^ : (x, n) < 0}) . 
So ([25]) shows {z4,n) ^ 0. Similarly, 6*14 < 624, implies {z^^n) < 0, a contradiction. Now 

cos ^23 — cos ^13 cos ^19 COS ^13(1 — COS ^12) tan(0i2/2) tan(l/20) 
cos B213 = ^ ^-7^ = — 7\ ^ 



sin0i3sin0i2 ^ sin0i3sin0i2 tan0i3 tan(7r/2 + 1/10) 

So, 9213 ^ arccos( ^Jg)^jffy\o) ) < Tr/2 + .01. By m, (^13 - 0i2? ^ Of^ ^ {Ou + Orsf, so 
7r/2 - 1/10 ^ 6*14 ^ 6*13 + 6*12 it/2 + 1/10 + 1/10. So we similarly conclude that G214 ^ 
arccos(4^^^) < vr/2 + .02. By m, 

cos 6*34 = cos(0i3 + ^'14) + sin 013 sin 014(1 + cos(27r - G213 - ©214))- (59) 

Since 6314 < vr by Proposition 13. 1 ( e) . 

|1 + cos 6314! = |l + cos(27r- 0213 -6214)1 ^ |l + cos(27r-7r-.03)| < .0005. 

Then ([58]) applied to ([59]) gives a dichotomy: (i) if 6*13 + 6*14 ^ vr, then 16*34 - (0i3 + 6*14) | < .04. Or 
(ii) if 013 + 6*14 > TT, then |034 - {2tt - 0i3 + 0i4)| < .04. Incase (i), 034 ^ 7r/2 + 7r/2 - 1/10 - .04 > 
TT - 1/2. In case (ii), 634 > 27r - {tt/2 + 1/10 + Tr/2 + 1/5) - .04 > vr - 1/2. So, in either case. 
Lemma 13.121 is violated. Therefore Case 1 cannot hold. 

(Case 2) 7r/2 + 1/10 < max{0i3,023} ^ 2.9. Without loss of generality, label the edges so that 
013 ^ 023, so it/2 + 1/10 ^ 013 ^ 2.9. By the triangle inequality on the sphere, 023 ^ 0i3 — 0i2 ^ 
013 - 1/10 ^ tt/2, so 

/ cos 012 - cos 013 cos 023 V „ /cos(l/10) - COS 0i3 COS 0i3 \ TT 

013 COS 6132 = 013 —7. r-7 ^013 —7. TTTTv^ > o- 

V sin0i3Sin023 / V sin0i3Sin(0i3 - 1/10) / 2 

Since 023,013 ^ tt/2, (l39|) gives 



^34 = ^Jel^ + 023 + 2013023 COS 0132 > ^2(^2)' + 2(7r/2) (^/2) = vr, 
so ()44p is violated. Therefore Case 2 cannot occur. 

(Case 3) The case max{0i3,023} > 2.9 cannot occur, by Also, max{0i3, 023, 0i4, 024} < tt/2 

implies that T4 and T3 are contained in separate but adjacent quarter spheres, so r4UT3 is contained 
in a half sphere. To see this, consider the edge Ci T3. Let ni2 be the plane such that ni2 Ci S'^ ^ 
T4 n T3. Temporarily suppose 0i3 = 023 = 7r/2, and let R: — )■ be the origin fixing rotation 
that maps the midpoint of T4 n T3 to V3. Then is contained in a quarter sphere whose boundary 
is contained in ni2 H 5^ and {R{Ili2)) H S^. If 0i3,023 ^ tt/2, then T4 remains in this quarter 
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sphere, since otherwise ?;3 ^ {x S 5^: dg2{x,vi) ^ 7r/2, ^52 (x, 112) ^ 7r/2}. So, Proposition I3.6r a) is 
violated. Ah cases therefore produce a contradiction, concluding the proof that 612 > 1/10. 

We now prove that sinBjjfc > 1/10. We argue by contradiction and assume that sin ©132 < 1/10. 
The procedure is similar to the calculations above. The theme is: if some Qijk is very small or very 
large, this creates additive relations among the edges, up to small errors. Before we begin, we list 
a few consequences of the spherical law of cosines, for i,j,k distinct elements of {1,2,3,4}. 

cos6'ij = cos(6'ifc -9kj) +sin6'ifcsin6'fcj(cos6ifcj - 1). (60) 

cos6'ij = cos(6'jfc + 0kj) +sin6'jfcsin6'fcj(cosGifcj + 1). (61) 

Now, assume 6132 < vr/2. Since 6132 + ©234 + ©134 = 27r and —©134 > — vr from Proposition 
Eme), ©132 + 0234 ^ vr. So ©234 > 7r/2 and sin ©234 ^ sin(7r-©i32) = sin(©i32) < 1/10. So, byre- 
labeling edges if necessary, we may assume that ©132 > 7r/2. Since ©132 > vr/2, cos ©132 ^ —.9949. 
Without loss of generality, label the edges so that ^23 = iiiax{0i3, ^23} and ^13 = min{0i3, ^23}- 
(Case 1') Assume max{ 6113, 6*23} < vr/S. By ([39]), ^'34 ^ 6*13 + 6*23 ^ 27r/3. But then Proposition 
13.6( a) is violated. (Let v G TiriT2 satisfy ds2{v,Vi) = min{034, 7r/2}, and note that {vi}f^i C {x G 
5^: d{x,v) ^ 7r/2}.) So Case 1' cannot occur. 

(Case 2') Assume max{6li3, 6^23} ^ vr/S and 1/10 < min{6li3, 6I23} ^ 1/2. Then 6*13 + 623 < n 
by Lemma [3.12( so ([571) and (f58]l give \6i2 — (^13 + ^23)! < -H- By the spherical law of sines 
and Lemma [3.12| sin ©123 = sin ©132 sin 6*13/ sin ^12 ^ sin©i32sin(l/2)/sin(l/2) < 1/10. Since 
^13 ^ ^23 and ^13 ^ 1/2, Lemma [3.121 gives cos 6*13 ^ |cos023| ^ cos ^23 cos 6*12. So, the formula 

cos ©123 = (cos 6*13 — cos 023 cos 0i2)/(sin 012 Sin^23) shows cos ©123 ^ 0, i.e. cos ©123 > .9949. 

Since 023 < — 1/2 from Lemma l3.12| 0i3 + 023 < vr. Combining this with cos ©132 < —.9949, 
(t6T]) and (jM]), we see that 0i2 ^ 013+023-.11. Since cos ©123 > .9949, ^ and 0i2 > .1 imply 024 > 
V(7r/3)2 + (.1)2 + 2(.l)(7r/3)(.9949) > tt/3. Now, since cos ©123 > .9949, ©123 + ©124 + ©324 = 27r, 
and —©123 > — vr from Proposition 13.1( e). ©123 + ©324 ^ vr, so cos ©324 ^ cos(7r — ©123) < —.9949. 
Since 023,024 > vr/3, by relabeling the edges if necessary we may assume that cos ©123 < —.9949, 
max{0i3,023} ^ 7r/3 and 1/2 < 7r/3 < min{0i3, 023}. Thus, we enter into the following case. 
(Case 3') Assume max{0i3,023} ^ tt/3 and 1/2 < min{0i3, 023}. Recall that we may label the 
edges so that 023 = max{0i3,023} and 0i3 = min{0i3, 023}. We claim that 023 > 1.3. If not, then 
(IHSD says 0I4 = 0^3 + 0^3 + 20230 13 cos ©132, SO 034 ^ Y^2(1.3)2 - 2(.994) (7r/3) (1/2) < 1.53 < 7r/2 
and {fi}f=i C {x G 5*2 : ds2{x,V3) ^ vr/2}. Thus, Proposition I3.6l fa) is violated, and our claim is 
proven. 

As in Case 1 in the proof that 9ij > 1/10, (j6ip gives a dichotomy 

(i) |012-(013 + ^23)| < .11, or (ii) |0i2 + 0i3 + 023 - 2vr| < .11 

In case (i), 0i2 ^ 023 + 0i3 — .11 > 1/2, and in case (ii). Lemma 13.121 savs 6ij < it — 1/2, 
so 012 ^ .89 > 1/2. So in either case, 1/2 ^ 0i2 ^ vr — 1/2, using Lemma 13.121 again. By 
the spherical law of sines, sin©i23 ^ sin ©132/ sin 0i2 ^ (1/10)(1/ sin(l/2)) < 21/100. Similarly, 
sin ©312 < 21/100. 

Assume (i) holds. If ©123 < vr/2, then cos ©123 > .977. By (i), 0i2 ^ 1.3 + 1/2 - .11 = 1.69. By 

§B,0h = ^1^2 + + 2^12023 cos ©123, so 

024 > V(l-69)2 + (7r/3)2 + 2(1.69)(7r/3)(.977) > 2.72 > ^ - ^• 

contradicting Lemma l3.12i If ©123 ^ vr/2 then cos ©123 < —.977. Then (j6ip gives a dichotomy: 
either (i)' |0i3 - (0i2 + 023)1 < -26 or (ii)' |0i3 + 0i2 + 023 - 2vr| < .26. If (i)' holds, then 0i3 ^ 
012 + 023-.26, but (i) says 012 > 613+623- -H- So 0i3 ^ 013 + 2023-.37, i.e. 023 < -185, contradicting 
that 023 > 7r/3. So (ii)' must hold. But then (i) implies 20i2 ^ 0i2 + 0i3 + 023 - -H ^ 27r - .37, so 
012 ^ vr — .185, contradicting Lemma 13.121 Therefore (i) does not hold. 
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Assume (ii) holds. Suppose 6123 < so cos 6123 > .977. Then (pO|) and ([58]) show that 
1^13 - 1^12 - 02311 < -26. If 012 ^ O23, then {Ou - (% + ^23)! < -26. So by {ii) 26^2 ^ 0i2 + ^is + 
^23 — -26 ^ 27r — .37, contradicting Lemma 13.121 If ^23 ^ G12, then the same argument shows 
2^23 ^ 27r — .37, contradicting Lemma 13.121 We may therefore assume that cos0i23 < —-977. If 
COSG312 > .977, we get the same contradiction, so we may assume cos 6312 < —.977. In summary, 
cos 6123, cos G312, cos 6132 < -.977. 

For i,j distinct, i,j G {1,2,3,4} let Eij = Ti Ci Tj. Let so that Oi'ji = maxj^is, ^23, 0i2}- 
Let Hi'j' be the plane containing Eiij' and the origin, and let n be the unit normal to Hi'j' such 
that {n, Z4) < 0. We claim that the edges -E'i'4, -Ej'4 satisfy 

Ei'A, Ef4 C {x G 5"^ : ^ (x, n) ^ .65}. (62) 

Also, the union of Ei/^, Ej/^, Ei/ji forms a noncontractible loop in this topological annulus. 
To prove Ei/4 C {x G : (x,n) ^ -65}, it suffices to show 

Qj'i'A > TT — sin~'^(.65), or ^4/4 < sin~^(.65). (63) 

To see this consequence, let x £ S'^ he contained in the great circle containing £^4/4 and let 
{i',j',k} = {1,2,3}. Since {n,Vi') = and {n,Zk) > 0, (x,n) is maximized when ds2{x,Vii) = 7r/2 
and {x,n) ^ 0. For such an x, (x,n) = sin0j/j/4. So, in the case that Qj'i'4 > tt — sin~^(.65), 
Ei'i C (x G S*^ : (x, n) ^ .65}. Let now x G {y G 5"^ : dg2{y, Vi') ^ 61/4} D Ei'4. For such an x, (x, n) 
is maximized for x in the plane containing {uj/,n, 0} with dg2{x,Vii) = Oin. So {x,n) = sin0j/4. 
So, in the case that 0j'4 < sin"^(.65), £"4/4 C {x G S"^ : (x,n) ^ -65}, as desired. The contain- 
ment Eji4^ C {x G S"^: (x,n) ^ .65} follows similarly if one shows: Qi'j'i > tt — sin~^(.65) or 
6j>4^ < sin~^(.65). Specifically, in the above paragraph we switch all indices of the form i' to j' and 
all j' to i', and the containment follows. 

We now discuss the proof of ([63]) . Let {i' ,j' , k} = {1, 2, 3}. Since Qj'i'k, &i'j'k < tt by Proposition 
13.1( e) and Ei'j/ C {x G 5^ : (x,n) = 0}, we conclude that Ei'k,Ej/k C {x G S"^ : (x,n) ^ 0}. By 
Lemma ETTll Y^j-, j-^i/ ex.p~^{j) = and Xlj: j^j/ exp~/(j) = 0, so we conclude Eii4,Eji4 C {x G 

: ^ (x,n)}. By the definition of Oi'j' and Lemma 13.141 Qj'i'^ ^ &ki'4- Ii^ particular, since 
Qj'i'A + 0fci'4 ^ we conclude that Qj'i'4 ^ vr/2. Similarly, Qi'j'i ^ -k/2. If £4/4 U Eiij' U Ej'4 were 
contractible in {x G 5^ : ^ (x,n) ^ .65}, then by (f63|) . the location of V4 is such that avoids 
the region {x G S'^ : {x,n) ^ .65 or (x,n) ^ 0}. But then either Qi'j'4 or Qj'i'4, is less than tt/2, 
a contradiction. So, the noncontractibility property of Ei'4 U E^'j' U Ej/^ follows from (j63p . which 
remains to be proven. 

Consider ^j/^ and Ei/ji as fixed with Qj/iik,Eii4 variable. From Lemma [3.14| recall that exp^^(4) 
and expT7^(/) have angle Qj'i'4. Also exp^^{k) and exp^^{f) have angle Qj'i'k, and expT7^(-) is 
confined to a plane. By Lemma 13.14^ scaling 9i'j' and 0j/fc by the same constant (greater than one) 
leaves Qj'i'4 unchanged and increases 6i'4. So in proving ()63p . we may assume that 6i'j' = tt — 1/2, 
by Lemma 13.121 

Suppose 9i'k < .815(7r — 1/2). In this case, we will show that sin0j/j/4 < .65, recalling that 
Qj'i'4 > 'tt/2. Using planar geometry and Lemma 13.141 one can show: decreasing Qj'i'k decreases 
Qj'i'4. So, for the purpose of showing sin0j/j/4 < .65, we may assume that cosGjVfc = —.977. 

Now, suppose 9iik ^ .815(7r — 1/2). In this case we show that sin0j/4 < .65, recalling that 
Oi'4 < vr/2, since ([I2]) says 

ef,^ = {e,,j, - Oi^kf + 2ei'j'ei'k{l + cosGj^fc) < (.185(7r - 1/2))^ + 2{tt - 1/2)2(.815)(.033) < .62. 
Consider 6i'k,Eirjr as fixed with ©j'j'fc, -E'i'4 variable. By Lemma 13.141 

^i'4 = I |exp,vH4)| I2 = I |exp^7^(/) + exp~^(A;)| = 9i,j, + ef,,, + 29i'j'ei'k cos Qj'i'k- 
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Recall cosQ j'i'k < —.977, so decreasing Qj'i'k increases 9i'4. So, for the purpose of maximizing 
sin 6*4/4, we may assume that cosQj'i'k = —.977. With these reductions (i.e. that cosGjVfc = —.977 
and 9i'j' = n — 1/2) for these two cases, one can then verify (j63p directly as a one-dimensional 
inequality, using Lemma 13.141 and treating Oi^ as a variable with all other quantities a function of 
6ik. With identical reductions, Qi'j'4 > vr — sin~^(.65) 9j'4 < sin~^(.65). Therefore, (|62|) holds. 

Given (|62p . note that the circumference of the geodesic ball {x G S'^ : {x,n) = .65} is given 
by 27rY^l - (.65)2 > 4 77^ gy ^he noncontractible property, 6*^/^ + 6*4/4 + 0j'4 ^ 4.77. By ([39]) . 

6*4/4 ^ {9i'j' - Oi'kY + '^(^i'j'(^i'ki'i- + cos Qj'i'k) ^ 6*4'/ - 6*4'^ and similarly 6*^/4 ^ (6*4/^/ - 9fk)- 
Substituting these inequalities into 9iiji + ^4/4 + 9j'4 ^ 4.77 gives 

9i'j' + {9i'j' — 9i'k) + {9i'j' — 9jik) ^ 4.77. 

However, by (zi), 9i'ji + 9i'k + 9j'k ^ 27r — .11. So adding the inequalities gives 404/j/ ^ 4.77+27r — .11, 
i.e. ^j/j/ > 2.73 > TT — 1/2, violating Lemma [3.121 So, all cases produce a contradiction, and the 
result is proven. □ 



4. Numerical Computations 

We now provide a more comprehensive analysis of the system of equations resulting from Lemma 
13. 8i We begin by writing this system explicitly. For i £ {1,2,3} choose j,i such that = 
{1,2,3}. Define 

rr , sin 6i£ s'mOji sin 9ij 

74 = - V A {vi,V4} = — — cos 6*4^ + + cos 6*4^. (64) 

Uie Vj£ Uij 

Then 7 = 72 from (j35p . Our system follows by cyclically permuting the 9ij in (\37\i . i.e., 

VXcos (^/9'f2 + 6*^3 + 26*126*13 cos 6213) +71 = 0, 
< \/A cos f V^12 + ^23 + 2012^23 COS 6123) + 72 = 0, (65) 
\/A cos ( V^1^3 + ^23 + 2^13^23 COS 6132) + 73 = 0, 

where Qije and A are defined in (jl2p and (|32p respectively. 

Observe that (f65|) gives three equations in the three unknowns 6*12,^13,^23. Define H: [0, vr]^ — )• 
so that i/(6*i2, ^13, ^23) = {Hi, H2, H^) is the vector corresponding to the entries of the left side 
of ([65]). We now wish to find (6*12,6*13,6*23) such that H = (0,0,0). 

We first examine neighborhoods of the two points discussed in Section [3j 

Lemma 4.1. The global maximum of F does not occur in the I2 halls of radius 1/100 around the 
points {6*12,6*13,6*23} = {9,9,9} £ ^^^f^ q ^ arccos(-l/3) and 6* = 1.53796841207904. 

Proof. Let 6 = 1/100. We claim that | (VksFq (6*12, 6*13, 6*23), ^ 20, for (6*12,6*13,6*23) in an £2 ball 
in of radius 5 ball around the point (arccos(— 1/3), arccos(— 1/3), arccos( — 1/3)) G M^, where 
V G 5^. Now, i<o(arccos(— 1/3), arccos(— 1/3), arccos (-1/3)) = 6(arccos(-l/3))2 ^ 21.9031 while 
F at the spherical propeller partition evaluates to (9/4)71^ ~ 22.206. So, given this claim, since 
(5|(V]R3Fo)| < 22.206-21.903 in 5|({6*}f^i, (5), no global maximum of F may occur in B|({6*}f^i, 5). 
Similarly, we claim that |(Vi[j3Fo(0i2, ^13, ^23)1 ^^)| ^ 15 for (6*12,6*13,^23) in an £2 ball of radius 5 
around {6*}^^^ £ with 6* ^ 1.53796841207904. Since Fo({6*}f^J « 21.7391, once again, no global 
maximum of F occurs in -B2({^}f=i, 5). 
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Proof of claim: for = arccos(— 1/3) or 6 = 1.537968, we first show that ||Hessian(i<o)||2 < 488 
for ah points in Bl{{9}f^^,6) C R^. Let F{x,y,z) = Fo{x,y,z). Then 

cos (y) y 



6 + 4. 



sin (y) 



cos (y) xy cos (x) (cos {z) — cos {x) cos (y)) xy (cos (2;))^ 

^ ^ ; ^ : ; ^ ~t~ 4 



sin {x) sin (y) 



^ (cos (z) — cos (x) cos (y)) y cos (x) ^ ^ 



(sin (x))'^ sin (y) 
(cos (2) — cos (x) cos (y)) xy 



cos (x) zy 



+ 4 



cos (z) z 
sin [z) 



(sin(x)) sin (y) 
„ cos (z) xz cos (x) 



sin (x) sin (y) sin (z) sin (y) 

(cos (y) — cos (x) cos [z)) xz (cos (x))^ 



sin (x) sin (z) 



xy 



(cos (y) — COS (x) COS (z)) z cos (x) 
(sin(x)) sin(2;) 

- 2xy o (^°^(^))^^^ I o^"'^^^^^ 



(sin {y) y 



sin (y) 



(sin (x))^ sin (z) 
(cos (y) — cos (x) cos (z)) xz 
sin (x) sin (z) 

^ (cos (x))^xy 
(sin (x))^ 



+ 2 



(cos (z) — COS (x) COS (y)) xy cos (x) cos (y) (cos (z) — cos (x) cos (y)) x cos (x) 



(sin (x))^ (sin {y)f 



(sin(x)) sin (y) 



cos (x) y (cos (z) — cos (x) cos (y)) y cos (y) cos (z) — cos (x) cos (y) 

~t~ 2 ; ; ; 2 ■ ; ■ ■ —75 -\- Z • 



+ 2 



sin (x) 
sin (x) zy cos (y) 



sin (x) (sin (y)) 

sin (x) z sin (y) xz cos (x) 

2 . , .... +2 



sin (x) sin (y) 
sin (y) z 



sin (z) (sin (y))^ sin (z) sin (y) (sin (x))^ sin (z) sin (x) sin (z) ' 

With 6 < 1/100, 9 = arccos(-l/3), and (x,y,z) € ({^}?=i> '^)> we have l/|sin(x)| < 1.1, and 
similarly for y, z. Therefore, in this 5-baIl we have \Fxx\ < 228. Here we have bounded the cosine 
terms by 1, the x,y,z terms by 2 (since arccos(l/3) + 6 < 2), and the inverted sine terms by 1.1. 
Using these same bounds, we find that \Fxy\ < 117 in this same (^-ball. For a 3 x 3 matrix A 



let 



=1 '^ij)^^^- Since F is symmetric in its arguments, we have the following bound in 



||Hessian(F)||2 ^ 



'228 117 117^ 
117 228 117 
117 117 228> 



< 



as desired. And the same estimate of the Hessian applies for 9 ~ 1.5379684120790425. 

Apply now the Hessian estimate to ||V]g3/(fe)||2 ^ \\'^R^fi'^)\\2^^^^Px£B'^{a,5) l|Hessian(/)(x)||2, 
which holds for every sufficiently smooth / : M'^ — t- M if a, 6 € satisfy ||6 — ^ 6. We then 
explicitly calculate | |Vk3Fo||j,j.c^os(-i/3)}3^J b < 13-6 and 1 1 VusFqI 115379684)3^ J b < 9.7. So, in 
S|({arccos(-l/3)}f^i,(5) we have IIV1R.3F0II2 < 13.6 + 4.88 and in S|({1.5379684}f^i , 5) we have 
IIVJ53F0II2 < 9.7 + 4.88, proving the claim. □ 

Lemma 4.2. We may assume that ^/X > .18. 

Proof. Substituting in ([H]), det(t;2, vi, V4) = det(t;i, ^2, ^3)^^^^. Using (j44l). 

sin(l/2) 



|det(?;2,^;i,'y4)| > |det(vi,'y2,'y3)| ^ . 

VA (vr - 1/2) 

Now, suppose y/X < sin(l/2)/(7r - 1/2), and \^ < sin(l/2)/(7r - 1/2), with A' defined by 



(66) 



A' 



sm &12 sm t'24 sm C714 
— V4 H vi H V2 



'12 



724 



'14 



23 



Applying ([66]) twice (with indices changed appropriately) shows |det(t;2, vi, ■^4)! > |det(t;i, t;2, fs)] 
and |det(t;2, f 1, ^4)1 < |det(t;i, ?;2) ^^3)!, a contradiction. We may therefore assume that one triangle 
of P (which, up to relabeling, is T^) satisfies ^/\ > (sin(l/2)/(7r - 1/2)) > .18. □ 



We now derive some modulus of continuity and derivative estimates for our system (I65p . Recall 
that P = {Ti}f^^ satisfies Proposition IS.lT e). Since (|65p must hold for the data of each individual 
triangle Tj, we select in ()65p the data from the triangle Then 612 + ^23 + ^13 ^ 27r, since the 
perimeter of T4 is bounded by 2tt. By relabeling the edges of T4 and using the pigeonhole principle, 
we may assume 612 ^ 27r/3. The following summary of ()65p holds for any labeling of the form 
= {1,2,3}. 



triangle inequality on the sphere 



6*12 + 6'23 + Oi3 < 27r 
012 ^ 2^/3 

^13, ^23 ^ TT- 1/2 

maxijg{i^2,3},i^j{^ij} ^ TT^-^ 

9^j > 1/10 

sin Qiji ^1/10 
> .18 



CoroUary[3lI] 
Lemma 13.161 
Lemma 13.161 
Lemma 14.21 



(67) 



For a map / : {X, 



ly) between two normed linear spaces, a modulus of continuity 
for /, denoted by dj, satisfies ||/(a;i) — /(x2)||y ^ C(Jj(||xi — X2\\x ,xi,X2)- We allow ujf to depend 
on xi, with X2 in a (5-ball around some given point xi = (^121^131^23)- 

Lemma 4.3. Suppose 6 < 1/100, and consider Hi: ([0, vr]^, 1 1-| |^) R, Hi = Hi(" " " ^ 
Suppose also that A = A(^i2, ^13, ^23) > V- Then, in an £qq ball of vadius 6 centered at ( 
the following holds 



'121 (7131 ^23;- 
^12i ^13i ^23); 



ujHAS)^G{6,e,r])=6 
+ {O12 + 6){6i3 + 6) 



7 15^ 

2 ^ 



155 



+ 3 2(^12 +5) + 2(^13 + 5) 



+ 



+ 



1 



^sm{9i2 + d) sm{9i3 + 6) sin(6'i2 + 5) sin((9i3 + 5) 
Or, applying Lemma \4-^ for 6 sufficiently small. 



ojhAS)^S 



7 



+ 



15 



2 2V.0329 - 155 



+ 3 2(^12 + 5) + 2(^13 + 5) 



+ (^12 + 5)(^i3 + 5) 



2 



+ 



+ 



sin(6'i2 + 6) sin(6li3 + 5) sin(6'i2 + 5) sin(6ii3 + 5) 



to 



Proof. We want to compute the modulus for the following function from 

^1(^121^131^23) = VAcos (^0^2 + ^13 + 2^12% cos 6213) +71 

where A and 71 are defined in (j33p and (j64p respectively. In what follows, we use the product rule, 
chain rule, and sum rule for the modulus ujf^.y 



(5)^a;^(5) + ||^/A|U.a;_^^ 



(6) + u;^AS) 



1 7 

^ ^Va('^) + 3 • ^ • ^{922+923+2012613 cos 0213) ('^) + o*^- 
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Here we used A ^ 9, \/A 3, |^(cos(\/t))| ^ 1/2, \{smt)/t\ ^ 1, and |^((smt)/i) 1/2|. Also, 
the term (7/2)5 appears since ^ 1 1 Vk37i(0i2, ^is, ^23)! Ii ^ 7/2- Now, observe that 



dX 



sin ^12 / d sinx 



89 



12 



712 



dx X 



sm Ui3 sint/23 ^ . 
+ 2 TT^^ (- sm6'i2j 



d sinx 



dx X 



x=9i2 / 
sin6>i3 

^13 



713^23 



cos 6^23 + 2 



d sinx 



dx x 



sm t/23 



723 



cos ^13. 



So, taking absolute values, and using that |sin(t)/i| ^ 1 and | Jj(sin(i)/t)| ^ 1/2, we see that 
\d\/d6ij\ ^ 5. So, if X{6) = A(0i2, ^13, ^23) > ?? > for some r/, and if x = {x,y,z) satisfies 
\x — 012 1 < 5,\y — 013 1 < 6,\z — 023 1 < then at x, 



5(^/A) 



90 



12 



1 1 



max 



Therefore, for such x. 



5A(0) ^ 1 

5012 2 y'rf 



1 



|Vk3A||^(^ 



1 

• 5 ^ - 



1 



2 V^T^TM 



• 5. 



I|Vm3(^/A)||i ^ 



15 



2V?? - 15c) ■ 

Combining this estimate with our estimate above, we get the following modulus estimate, that 
applies to x, as just described 

15 



7 



2 2Vr/ - 155 
+ {912 + S){ei3 + S) 



+ 3( 2(012 + 5) + 2(013 + 5) 
2 



+ 



+ 



{91 
1 d^p 



2 89 



sin(0i2 + 5) sin(0i3 + 5) sin(0i2 + (5) sin(0i3 + 5) 
! function x 1— t- x/sinx a 

(cos 023 - cos 012 COS 013)), 

'sin 012 COS 013) 



We used here the monotonicity of the function x i— t- x/sinx and the following identities for 

m2,^13,^23) = (e?2+^?3 + 2 



Sin &12 sin &13 
^12^13 



12 



^12 + - 

sinWi2 sinWi3 

+ (cos 023 - COS 012 cos 0i3)(0i2/ sin 012)' • (013/ sin 013) 



3 , Pl2fl3 a , r> / 1 
712 + ■ ^ COS 013 + COS 9213 1 



1 di) 
2d9~ 



h3 + 



sinwi3 
'12^13 



712 



tan 012 



713, 



(Sint/i3 COSt/i2j 



13 



sinfi2 sinWi3 

+ (cos 023 - COS 012 COS 013)(013/ sin 013)' 



712 



/ sin 012) 



1 di/j 
2d9~ 



23 



a , '712(713 a , r> / 1 
c'13 + ■ o ^12 + COS 6213 1 

sin0i2 

^12^13 ^ 
sin 012 sin 0i3 



713 



tan 013 



712, 



sm 6/23 J 



□ 



Analogous estimates give a gradient bound for Fq, as defined in (j43p . Since Fq is symmetric in 
its arguments, it suffices to take a derivative with respect to one variable. Fix (012,^13,^23), and 
let < 5 < 1/100. For any (x, y, z) such that |x — 0i2| < 5, |y — 0i3| < 5, |z — 023 1 < 5 we have the 
following derivative bound 



25 



9Fn 



de 



12 



. 2(% + ^)(g23 + <?) , , (gl2+<^)(% + <^) 

" ^ ^ sm(^i3 + 5) sm(023 +5) sin(0i3 + 5) 



(70) 



sin(023 + 5) ' 'V tan(0i2 + (5). 

The restrictions of (j67p and the bounds of (j68p and (|70p finahy yield Theorem II. li 

Proof of Theorem \l.l\ Let P be a partition of the sphere S'^ that maximizes Ylt=i ll-^jlli- Arguing 
by contradiction, we may assume that this partition exceeds the value of the spherical propeller 
partition, 97r^/4. From Proposition 13.11 we may assume that P is a collection of four nonempty 
spherical triangles along with four vertices such that det(f i, 112, 't's) > 0. Since P is a global 
maximum of P, it satisfies three equations in the three variables {^12, ^13, ^23} defined by (1650 . We 
may restrict these variables via (I67p . We may also exclude neighborhoods of the two candidate 
maxima found in Section [3l using Lemma 14. 1[ Combining the estimates of (I68p , (I69p and ()70p gives 
sufficient information for an e-net traversal of the parameter space defined by (I67p . This numerical 
computation finds that no zeros of the system (1650 . within the parameter space defined by (I67p . 
exceed the objective value of the propeller partition. This computation contradicts the assumed 
existence of a maximal P with four nonempty elements. So, recalling ([8]), X^f^i \ \zi\\2 ^ 97r^/4 = 
27r3(9/(87r)). □ 

4.1. Implementation. We describe the e-net traversal of Theorem ll.ll further below. Before doing 
so, we describe some issues related to the numerical verification of inequalities. In particular, we 
discuss how a computer can rigorously verify a mathematical statement. Recall that a normal 
double precision floating point number in the IEEE- 754-2008 standard is a number in base two of 
the form 

±l.aia2a3---a52 x 2^1^2-611-2"'^ 

where Oj E {0, 1} and 61&2 • • • &11 G [2, 2^"^ — 1] n Z. All of our computations are done with this 
discrete set of numbers, with zero included as well. Note that the distance between two consecutive 
numbers changes according to the values of the two numbers. For example, the distance between 1 
and the next largest number is 2~^^, whereas the distance between 2 and the next largest number is 
2~^^. (In general, if the number 2^' and the next largest number can be represented in this system, 
then their distance is 2~^^~^^.) Due to this spacing between numbers, one must round the result 
of any arithmetic operation. For example, the addition 1 -|- 2~^^ evaluates to 1, if we choose to 
round to the nearest double precision number. For the sake of flexibility, our analysis below applies 
regardless of the rounding scheme that is chosen. 

For x,y (z [0,2^*^^] where x,y are normal double precision numbers, let fl(x -|- y) denote the 
normal double precision number that is evaluated in a computer with the IEEE-754-2008 standard. 
It suffices for now to note that {x + y) ^ {1 + 2 ■ 2~^^)fl(x + y). In this case, we say that {x + y) 
has multiplicative error bounded by (1 -|- 2 • 2^^^). When we speak of numerical error, we mean 
multiplicative error, unless otherwise stated. Define 

£ = 3-2-52^ 

To verify inequalities in our numerical computations, we perform operations for several floating 
point numbers and keep track of error bounds of the form [x + y) ^ (1 -|- 2 • 2~^^)fl(2; -|- y). Since 
(1 + 2 • 2-52)7V ^ (1 + AT . 3 . 2-52) = (1 + iVe) for 1 iV ^ 30000, we can bound the multiphcative 
error of N operations with a multiplicative term of the form (1 -|- Ne). (Here and below we are 
overly conservative in our error estimates.) 

We now present the formulas that are used in our computations, delaying further discussion 
of numerical errors. Let x,y > be two normal double precision floating point numbers with 
x > y. A loss of significance refers to a loss of binary significant digits in the computation of 
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fl(x — y). One can show that fl(x — y) loses at most q and at least p significant binary digits, where 
2"'? ^ 1 — (y/x) ^ 2~P. Due to potential loss of significance errors, we evaluate A as 



A 



^23 



sin^ ^13 ^ sin^ 6*12 



92 



'12 



+ 2 



sin f 23 sin fc/i3 



723^13 



cos 012 



sm «23 sm t/12 „ , sm t/12 sm t/i3 
H 7r~^ cos 6113 H — cos 623 



(71) 



P23'712 C712'713 

Our numerical computation of inequalities consists of four steps. 
(Step (i)) Suppose we want to check ([67|l for a closed ioo ball (i.e. a cube) of radius S with center 
,0i3>^'23) G [0, vr]'^. To do so, we must verify the following perturbation of (i67|) . 



n2, 



(% + Oje + 3<5) ^ 

Ou + 023 + ^13 ^ (2vr + 35) 

012 5^ (2^/3 + 5)(l + 10e) 

013 + 1/2, 023 + 1/2 ^{tt + 6){1 + We) 
{5 + max,jg{i^2,3},i^j{^ij})(l + 7e) ^ 
(0,j+5)(l + 7e) ^ (l/10)(l-e) 



(1 - 5e) 



(35 + ^(sin0ij sm9ji + 26)){1 + 1500e) ^ (cos0i^ - cos0ij cosOje + 2000e) 



(72) 



(cos ( 



cos t>ij COS Uji 



2000e) ^ (35 + ^{sm9ij sinOji + 25))(1 + 1500e) 



(155 + 2 • lO^e + X)> .03293 

If a given inequality of (j72p does not hold, then the corresponding inequality of (j67p is violated for 
the entire £00 ball of radius 6. 

(Step (ii)) If all of the inequalities of ([72]) are satisfied, we then check if (0i2, 0i3, 023) lies in a 
neighborhood of two of our known candidate global optima. Let x = {xi,X2,X3) G [0,7r]^. Suppose 
xi = X2 = X3 = arccos(— 1/3) or xi = X2 = x^ = 1.5379684120790425. From Lemma l4.ll we know 
that if = (012, 013, 023) is close to x, then F is bounded by (9/4)7r2 in !b^((0i2, 0i3, 023), 5). That 
is, if the following is not satisfied, then F is bounded as just stated. 



(d^2 (0, x)y + 40000e + 25 df,^ (0, x) + 6000e + 35 (1 + 20e) > ((1/100)2)(1 - 2e) 



(73) 



Here di^ and di-^ denote the usual £2 and ii metrics on M.^. 

(Step (iii)) If the verification of (|73|) fails, we use the modulus estimates of Lemma [4.31 to see if H 
is near zero at 0, with H = {Hi, H2, H^) defined after (f65|) . Let a, 13, a G {±1}- We first check if 
the following expression has constant sign, over the 2^ such choices oi a, (3, p,a. 



cos W0f2 + 0^3 + 2012013 



COS 023 — COS 012 COS 013 + 200060; ' 



sm ^12 sm ^13 



+ 25000e/3 



(74) 



\/A - 2 • 10%o- + (71 + 3000ep). 



If (|74p is well-defined and its sign is constant, we then consider the following inequalities involving 
H for i = 1 and (168)1. 



mm 



COS \ 01^ + 61^ + 26126 



13 



COS 023 — COS 012 COS 013 + 2000ea \ 



+ 25000e/3 



VA - 2 • lO^ecj + (71 + 3000e/9) 



sin 012 sin 013 
(1 - 3000e) > (1 + 3000e)G(5, 0"*, A - 2 • lO'^e). 



(75) 
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mm 

a,(3,p,ae{±l} 



COS I J 9^^ + 91^ + 2012 



'13 



cos 023 — cos 012 COS ^13 + 2000ea \ 
sin 6*12 sin ^13 J 



+ 25000e/3 



• \/A - 2 • 10%cj + (71 + SOOOep) 

7 



> (1 + 5e)V15(^ + 5(1 + 3000e) 

2 



(1 - 3000e) 

+ 3(2(^12 + <5) + 2(% + <5) 



(76) 



{9i2 + S){9n + 6) 



+ 



1 



mn{9i2 + S) sm{9is + 5) sin(6li2 + 5) sin(6'i3 + 5) , 
The inequalities for H2 and are constructed by permuting cyclically the indices appearing 
above. If at least one of (j75p and ()76p is satisfied (for at least one Hi), then Hi is nonzero on 

^((012, 013, 023), <5). 

(Step (iv)) If (|75|) and (|76|) are not satisfied for {Hi}f^^, we finally check the value of Fq directly, 
using a modification of (I43p . Define 



cos efj-fc 



^0^ 



Gf2(0,<5) =6(012 + <5) + 
+ 2 



cos 9ik — cos 0jj cos 0jfc + 2000e \ 
sin 6ij sin 0^,^. / 

3 {OI2 + 0^3 + Oh) + 2cos(e|i3)0l2013 + 2cos(ef23)^12023 + 2 COs(e|3i)023013, 
4(013 + <5)(023 +5) , J012+5)(013 + 5) 



sin(0i3 + J) sin(023 + 6) 

12 + 5)(023 + 5) 



+ 2- 



sin(0i3 + 5) 



sin(023 + 6) 



+ 2(013 + ^23 + 26) 



1 



+ 5 



(tan(0i2 + <5) 



;i + 1200e) + e 



If the following is satisfied, then Fo is bounded by (9/4)7r2 in B^{{9i2, ^13,^23),^), and therefore 
F satisfies the same bound for any global maximum of F in this ball. 

(9/4)7r2(l - 5e) - F^ • (1 + 5000e) > (5(G^2 + + G|3)(l + W^e) + lOOe. (77) 

Given these estimates, we can now describe the e-net traversal of Theorem 11.11 We need a 
discrete subset Af = {xi} C [0,7r]'^ together with a set of radii 6i > such that 

U Blix,,5i)^[0,7Tf, 

and such that one of the following cases occurs, for {xi,6i) with Xi G M. 

(i) One inequality of ([72]) does not hold. 

(ii) ([73]) does not hold. 

(iii) (f75|) or ([76|) holds, and (j7i|) has constant sign over all sign choices a,/3,cr,p G {±1}. 

(iv) dZZl) holds. 

In case (i), B'^{xi, 6i) is either outside of the domain of ([65]) or ([65|) cannot be satisfied in B^{xi, 6i). 
Case (ii) implies that F = J2i=i \\^i\\2 bounded by (9/4)7r^ on B'^{xi,6i). Cases (iii) and (iv) 
imply that no global maximum of F occurs in B^{xi,5i). So, the existence of M as defined above 
completes the proof of Theorem II. li 

To construct M, we perform the following standard procedure. Begin with the grid 

AAo = ^([0,100]3nZ3) C [0,7r]3. 

For each point x G A/q, check the inequalities discussed above with 6 = = 7r/200. If we discover 
that X satisfies one of (i) through (iv), then eliminate x from the set A/q. Next, include the element 
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(x, 6) in the set A^. Let A/'q be the set of all points x G A/q that have not been included in A^. We 
now "refine" the grid A/'q by a factor of 10 and repeat the above procedure. That is, we define 

M = (Y^([o,iooo]3nz3))n| J bE{x,,7t/200)] . 



For each point x G A/i, we now check (i)-(iv) with 5 = 5i = 7r/2000, and so on. For j ^ 1, at the 
j^^ iteration of this algorithm, we define 







X x^ 


(2l!)sin(x) 




sin(x) (21!) 



^3 = (Y^([0'io'^']'nz=')) n (J i?^(x„vr/(2 • lo^+i)) 

We continue in this manner until the procedure terminates. Upon termination at the j^^ step, 
= 0, so AA satisfies our desired properties and proves Theorem ll.il 

We now return to our discussion of numerical error. To keep track of rounding errors resulting 
from subtractions, we use absolute error where appropriate and add the e terms that correspond to 
absolute numerical errors. These additive terms are scaled according to the absolute value of the 
numbers being subtracted (due to the unequal spacing of fioating point numbers described above). 

In the description and analysis of (i) through (iv), 5 terms come from modulus of continuity 
estimates, and e terms account for numerical error. The latter is discussed further below. The 
first two inequalities of (I72p contain no e terms, since they can be checked in exact arithmetic. Let 
rr G [0, vr — 1/2 + 2 X 10"^]. (We make this restriction due to the fourth inequality of ([72]) .) Let us 
discuss the numerical computation of the sine function. If x > vr/2, we replace x with vr — x, and 
we evaluate sin(7r — x). Thus, we only need to consider x G [0,7r/2]. For such x, we first use a 19*'* 
degree (ten term) Taylor expansion around zero, summing terms in ascending degree. In perfect 
arithmetic, we get a multiplicative error of one plus 

Ti (7r/2)20 e 
s: ' ' < 2.6 X 10"^^ < -. 
2 21! 2 

Since x G [0,7r/2], there are no loss of significance errors in the summation of the series. The 
addition of each term of the Taylor series involves one operation, and there are ten such terms. 
Moreover, the n^^ term contains n + 1 operations. In total, one evaluation therefore has less than 
210 operations. So, very conservatively, the multiplicative error of sin(x) is bounded by (1 + 250e), 
if x is represented exactly. If the argument x has multiplicative error (1 + ke), then this same 
analysis shows that the multiplicative error of sin(2;) is bounded by (1 + (250 + 100A:)e), since the 
term x appears 100 times in the formula for sin. For example, since 9ij has multiplicative error 
(1 + 3e), we conservatively bound the error of sin^jj by (1 + 550e). 

To compute cos(x) for x G [0, vr — 1/2 + 2 x 10"^], we instead compute sin(7r/2 — x) with the 
procedure outlined above. Due to our representation of the points we compute tt/2 — x as a 
rational multiple of vr, so no loss of significance errors occurs in this subtraction. Hence, the error 
estimates for cos 6ij are exactly the ones already used for sin0,j. However, when evaluating the 
cosine of more complicated expressions, we need to revert to absolute errors. 

We will also need to use the square root function. Here we use the usual Newton's Method 
iteration, which is well known to have multiplicative error bounded by (1 + e). 

We now describe the error terms appearing in ()72p . which are perturbations of ()67p . The guiding 
principle is that we want to perturb each inequality so that it is more easily satisfied. To see why this 
is done, suppose an inequality of ()72p is not satisfied. Then, a modulus of continuity estimate and 
a numerical error estimate, used in the definition of (I72p . show that the corresponding inequality 
of ([671) is not satisfied in an £qq ball of radius 6. For example, contrast tlie first inequality of (I67p 
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with that of (j72p . We need only add an additional 6 term for each 9ij term, since the modulus 
of continuity of the function (x,y,z) i— t- x -\- y -\- z, viewed as a mapping from (M^, — M, is 

exactly 35. 

So, if the inequality of (j72p is violated at a given point 9 = (^125^13)^23)) then the modulus of 
continuity bound shows that the corresponding inequality of (|67p is violated in the ball B^(9,6). 
The same considerations apply to the second inequality of (I72p . Note that the first and second 
inequality are checked in exact arithmetic in our computer implementation, so we do not need to 
take into account numerical error. As we have mentioned above, our 9ij and 6 are represented as 
numbers of the form n^a/c) where a and c are integers, and it suffices to perform additions of the 
form 7r((a + b)/c), and the integer addition a + b can be performed with no error. 

The third inequality of (I72p only requires a modulus of continuity term of 5. However, we must 
now take into account numerical errors. We count ten numerical operations on positive numbers 
(including the error in representing vr as a floating point number three times), so in the worst case, 
the multiplicative error is bounded by (1 + lOe). Similar considerations apply in the next three 
inequalities. In the seventh and eighth inequalities, we use our analysis of the sine function, and 
recall that the multiplicative error for these operations is always bounded by 1 + 550e. So, in total, 
the evaluation of the term —(3(5 + ^^{sin 9ij sin 9 j£ + 25)) has multiplicative error bounded by 
1 + 1500e. Here we have also included 5 terms corresponding to the moduli of continuity, as before. 

Unfortunately, we cannot use multiplicative error to deal with the cosine terms, since loss of 
significance can occur in the evaluation of (cos 9i£ — cos 9ij cos 9je). However, since the cosine terms 
are bounded by 1 in absolute value, the absolute error from each cosine evaluation is bounded 
by 550e. Therefore the total absolute error from the right hand side of the seventh inequality is 
bounded by 2000e. This bound applies similarly to the penultimate inequality of ()72p . The final 
inequality of (j72p is similar to the previous cases. To get the 5 terms, we examine (|7ip and use 
that |(sin(a;)/a;)'| < 1/2 and |sin(x)/x| ^ 1 for x S [0, vr]. To get the e terms, note that each term 
s\v?9ij/9fj has absolute error bounded by 1200e. And each term sin 9ij sin 9 jk cos 9ik/9ij9jk has 
absolute error bounded by 1700e. The total absolute error of A is therefore bounded by 2 • lO^e. 

We now treat (j73p . First, recall that (jlSp yielded our two known zeros of the system (j65p . The 
first zero 9ij = 9 = arccos(— 1/3) is known exactly (up to multiplicative error (lie)). However, we 
can only find the second zero computationally. For the purposes of our computations, we therefore 
show that this second zero is contained in an interval of the form 9ij = 9 £ [x — 2000e, x + 2000e] 
where x = 1.5379684120790425. To see this, we take (jlSp and write it as a quotient of the left side 
over the right side, and we wish to find where the quotient is 1. For 9 € [x — 2000e, x + 2000e], 
the multiplicative error involved in calculating this quotient is bounded by 1 + 3000e. This follows 
by our analysis of the sine and cosine functions, and since this particular region of parameters 
avoids loss of significance errors. We can then check numerically that, for 9 = x — 2000e, this 
quotient, multiplied by (1 + 3000e) is less than 1, while for = x + 2000e, the quotient, multiplied 
by (1 — 3000e), is greater than 1. So, by our numerical analysis and continuity of the given function, 
a solution of (fTSj) must lie in the interval 9 £ [x — 2000e, x + 2000e]. 

We can now confront the errors of (I73p . Once again, the guiding principle is to add error 
terms so that the inequality is more easily satisfied. First, the term di^{9,x)'^ is the sum of three 
squared differences. Due to numerical error in the operations, and the uncertainty of the zero 9ij ~ 
1.5379684120790425, each squared difference has an absolute error bounded by 27r(2000e) + lOOe. 
(The factor of 27r comes from the squaring of the difference, and since ^ vr.) So, in total, the 
absolute error from the term d^^ {9, x)'^ is bounded by 40000e. Since we want to check this inequality 
in an £00 ball of radius 5, we use the modulus of continuity of d'j^ , which is 25d£-^ . We therefore add 

the term 25di^{9,x), and again add terms to this quantity to take into account its own absolute 
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error and its own modulus of continuity. The combined effect of summing the error terms then 
accumulates a multiplicative error bounded by 1 + 20e. 

We now explain the error terms of (j75p . From the analysis of the cosine, we know that cos 9ij 
for 6ij G [0,7r — 1/2 + 2 X 10^^] has an absolute error bounded by 550e. (This follows since the 
multiplicative error is bounded by (1 + 550e), and |cos(-)| ^ 1.) So, the expression cos023 — 
cos ^12 cos has absolute error bounded by 2000e. The remaining operations under the square 
root of ()75p all have multiplicative error bounded by (1 + 1200e). 

When the cosine of the square root term is evaluated, we need to again take into account absolute 
errors (since, unlike in the cos 6ij terms, we must now account for loss of significance). Since the 
argument of the cosine term is bounded by vr, we can bound the absolute error of the square root 
term by 4000e. To take the cosine of the square root term, we need to return to our Taylor series 
analysis. Note that our Taylor series has Lipschitz constant bounded by 6. (This is a conservative 
bound that does not account for cancelations in the series summation.) Therefore, the total absolute 
error of the cosine of the square root in (j75p is bounded by 25000e. From before, the absolute error 
in the computation of A in ()7ip is bounded by 2- lO^e. Similarly, the computation of 71 has absolute 
error bounded by 3000e. 

Temporarily ignoring any numerical errors, note that we are computing a bound for the absolute 
value of the derivative of a function. So, in order to conclude anything from (f75|) . the sign of (j7i|) 
must be constant over all possible numerical errors. (If not, it may be that (I74p has an actual value 
of zero.) Given that this sign is constant, we then take the minimum over all possible numerical 
errors. It then suffices to take the minimum over the extreme points of a suitable rectangle, due 
to monotonicity of (I75p with respect to the error terms. Note that we take the minimum since 
we wish to make the inequality more difficult to be satisfied. Finally, the addition of all terms on 
the left side of the inequality incurs a multiplicative error bounded by 1 + 3000e, and the addition 
of the terms on the right side has multiplicative error bounded by 1 + 3000e. A similar analysis 
applies to (f76]l . 

We now briefly mention the analysis of (|77p . which is routine at this stage. The principle that 
we apply here is to make (|77p more difficult to be satisfied, so we make its right side larger and 
its left side smaller. The multiplicative error of the Fq term is bounded by 1 + 5000e, since each 
cos Qijk term has multiplicative error bounded by 1 + 1200e. Each C-j term has multiplicative 
error bounded by 1 + 2500e, so the right side of (I77p has multiplicative error bounded by 1 + lO^e. 
A lOOe term is added to the right side of (I77p to account for the error of the subtraction on the 
left side. The 1 + 1200e term and the e term in the definition of G?- account for the errors of the 
tangent function, and the loss of significance in the subtraction that occurs in the square brackets. 

5. Discrete Harmonic Maps into the Sphere 

We now summarize and reinterpret the results of Section [3l Recall that we defined G{P) = 
J2i^i<j<i4^ij the sum of the squared lengths of the edges of our partition P. We then saw from 
Propositions 13.131 and 13.11 that all maxima of F are also critical points of G. We also found that 
critical points of G satisfy (j5ip . Given a partition P, think of its edges as rubber bands. (Recall 
that a rubber band has energy given by its length squared, and it exerts a force proportional to 
its length.) The critical points of G are given by configurations of rubber bands (arranged like 
the edges of Figure [2]^b) ) that are at equilibrium with respect to the forces exerted by the rubber 
bands. With this interpretation, (jSip becomes intuitively clear. This equation says that, at a fixed 
vertex, the (tangential) forces of the incident rubber bands are in equilibrium. (For the Euclidean 
version of this, see e.g. [10] and [3T].) 

These considerations lead us naturally to the following generalization. Let (V, E) be a (finite) 
graph. Let u: iy,E) — t- S"^ be an embedding of iV^E) into S"^ {u is sometimes called a map). To 
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be consistent with the above, let u{V) = {vi}, and label Length(«(ujf j)) = Oij where ViVj denotes 
a spherical geodesic connecting Vi to vj. Then a harmonic map u is a critical point of the energy 
functional E{i,j}&E^Ij- 
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